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Abstract 


The  recent  interest  in  the  development  of  small  UAVs  and  micro  air  vehicles  has  revealed 
a  need  for  a  more  thorough  understanding  of  the  aerodynamics  of  small  airplanes  flying  at  low 
speeds.  In  response  to  this  need,  the  present  work  provides  a  comprehensive  study  of  the  lift, 
drag,  and  pitching  moment  characteristics  of  wings  of  low  aspect  ratio  operating  at  low  Reynolds 
numbers.  Wind  tunnel  tests  of  wings  with  aspect  ratios  between  0.5  and  2.0  and  with  four 
distinct  wing  planforms  have  been  conducted  at  chord-Reynolds  numbers  in  the  range  of  70,000 
to  200,000.  In  addition,  the  effect  of  leading  edge  shape  and  fuselage  bodies  has  been  studied. 

As  an  example  of  an  application  of  this  wind  tunnel  data,  the  experimental  results  are  used  as 
part  of  an  aerodynamic  analysis  procedure.  This  procedure  is  incorporated  into  a  genetic  algorithm 
optimization  program  that  generates  optimum  MAV  configurations  given  certain  requirements  and 
constraints.  Results  obtained  by  use  of  this  optimization  procedure  have  revealed  that  useful  and 
accurate  design-optimization  tools  can  be  developed  based  on  the  experimental  data  presented 
within  this  report. 
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Chapter  1 
Introduction 


1.1  Motivation  for  Research:  Small  UAVs  and  MAVs 

A  great  deal  of  interest  has  emerged  during  the  last  half-decade  for  a  totally  new  kind  of 
unmanned  aerial  vehicle.  Although  large  unmanned  aircraft  such  as  the  Hunter,  Pioneer,  and 
Predator  have  been  successful  for  many  years  as  theater-level  military  reconnaissance  platforms, 
a  need  for  a  smaller,  platoon-level  unmanned  aircraft  exists.  Small  UAVs  such  as  the  Pointer, 
Sender,  and  more  recently,  the  DragonEye,  have  arisen  to  fulfill  this  need.  An  even  smaller 
alternative  is  the  micro  air  vehicle  (MAV).  MAVs  are  expected  to  be  carried  by  soldiers  in  the 
field  and  to  be  deployed  for  the  purpose  of  answering  the  question:  “what’s  over  the  next  hill?”. 
For  this  puipose,  MAVs  are  envisioned  to  be  small,  semi-autonomous  airplanes  with  limited  but 
versatile  performance  characteristics.  One  proposed  set  of  requirements  [1]  suggests  that  MAVs 
should  have  a  maximum  dimension  of  approximately  15  cm  (5.9  in),  be  easily  controlled  by 
soldiers,  and  be  able  to  operate  at  short  to  medium  ranges  in  the  order  of  10  Km  (6.2  mi),  with 
an  endurance  of  approximately  30  minutes.  In  addition,  MAVs  should  be  capable  of  carrying  a 
variety  of  sensing  and/or  surveillance  payloads  with  a  typical  payload  of  18  g  (0.63  oz).  Successful 
operational  MAVs  meeting  some  of  these  proposed  requirements  include  AeroVironment’s  Black 
Widow,  BAE/Lockheed  Martin’s  MicroSTAR,  and  the  Naval  Research  Laboratory’s  Mite  and 
Mite  n. 

The  design  of  efficient  MAVs  is  hindered,  however,  by  the  lack  of  a  thorough  understanding 
of  the  aerodynamics  associated  with  small  airplanes  flying  at  low  speeds.  The  operating  cruise 
speed  of  MAVs  is  typically  between  16  and  80  KPH  (10  and  50  MPH),  yielding  an  operating 
Reynolds  numbers  range  between  roughly  70,000  and  200,000.  Furthermore,  the  requirement  of 
smallness  and  compactness  dictates  the  use  of  very  low  aspect  ratio  (LAR)  wings  (approximately 
AR  =1).  A  need  exists  for  detailed  aerodynamic  analysis  tools  that  are  applicable  to  the  Reynolds 
number  and  aspect  ratio  operating  conditions  found  in  MAVs  and  in  some  small  UAVs. 

The  present  work  provides  comprehensive  aerodynamic  data  for  low  aspect  ratio  wings 
operating  at  low  Reynolds  numbers.  The  data  is  based  on  wind  tunnel  experiments  of  wing  models 
with  aspect  ratios  between  0.5  and  2.0  at  Reynolds  numbers  between  70,000  and  200,000.  The 
effect  of  wing  planform  is  also  studied  as  part  of  this  project.  Four  distinct  planforms  are  used  for 
all  aspect  ratios  and  Reynolds  numbers.  In  addition  to  this  basic  set  of  data,  experimental  results 
are  also  included  that  provide  general  information  about  other  aerodynamic  effects,  including 
the  increase  of  drag  due  to  a  fuselage  body  and  the  influence  of  leading  edge  shape  on  lift  and 
pitching  moment. 
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1.2  A  Short  Discussion  of  LAR  Wings  at  Low  Re 

The  aerodynamics  of  low  aspect  ratio  wings  at  low  Reynolds  numbers  have  essentially  not 
been  studied.  Investigations  of  the  two-dimensional  aerodynamic  characteristics  of  airfoils  at 
low  Reynolds  number  have  been  conducted  at  the  University  of  Notre  Dame  by  Pelletier  and 
Mueller  [2],  Three-dimensional  LAR  wings  have  been  extensively  researched  at  higher  Reynolds 
numbers  in  the  form  of  delta  wings  at  subsonic,  transonic,  and  supersonic  speeds.  Many  of 
these  studies  have  focused  on  the  high-angle-of-attack  aerodynamics  of  delta  and  other  types 
of  pointed  LAR  wings.  For  MAVs,  delta  wing  planforms  are  not  ideal  because  for  a  given 
maximum  dimension,  delta  wings  offer  less  lifting  area  (and  therefore  less  payload  capability) 
than  a  rectangular  or  circular  wing  of  the  same  maximum  dimension. 

Some  information  is  available,  however,  regarding  non-delta  LAR  wings,  with  much  of  the 
research  having  been  done  between  1930  and  1950.  Zimmerman  [3],  Bartlett  and  Vidal  [4], 
and  Wadlin  [5]  performed  experiments  of  LAR  wings,  although  at  significantly  greater  Reynolds 
numbers.  Theoretical  and  analytical  treatises  of  LAR  wing  aerodynamics  were  considered  by 
Bollay  [6],  Weinig  [7],  Bera  and  Suresh  [8],  Polhamus  [9,  10],  and  Rajan  and  Shashidhar  [11]. 

Perhaps  the  most  complete  analysis  and  review  of  LAR  wings  was  performed  by  Hoemer 
in  his  two-volume  series  on  lift  and  drag  [12,  13].  Hoemer  reviewed  many  of  the  theories 
developed  for  LAR  wings  of  non-delta  planforms.  A  variety  of  correlations  as  well  as  analytical 
methods  were  presented  and  compared  with  the  available  experimental  data  of  the  time.  Although 
the  information  presented  by  Hoemer  corresponds  to  higher  Reynolds  numbers  than  the  ones 
considered  in  the  present  work,  the  aerodynamic  theory  still  holds. 

This  theory  correctly  predicts  that  as  a  finite  wing  of  a  given  aspect  ratio  generates  lift, 
counter-rotating  vortical  structures  form  near  the  wingtips.  These  vortices  strengthen  as  the  angle 
of  attack  increases.  For  a  low  aspect  ratio  wing,  the  tip  vortices  may  be  present  over  most  of  the 
wing  area  and  therefore  exert  great  influence  over  its  aerodynamic  characteristics.  Wings  of  AR 
less  than  or  equal  to  1  can  be  considered  to  have  two  sources  of  lift:  linear  and  nonlinear  [13]. 
The  linear  lift  is  created  by  circulation  around  the  airfoil  and  is  what  is  typically  thought  of  as  lift 
in  higher  AR  wings  (Figure  1.1).  The  nonlinear  lift  is  created  by  the  formation  of  low  pressure 
cells  on  the  wing  s  top  surface  by  the  tip  vortices,  such  as  that  observed  in  delta  wings  at  high 
angles  of  attack  [14].  This  nonlinear  effect  increases  the  lift-curve  slope  as  a  increases  and  it  is 
considered  to  be  responsible  for  the  high  value  of  ostall  (Figure  1.2). 
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Pure  circulation 


Incompressible  flow  over  an  airfoil 


Figure  1.1:  Linear  Lift  of  LAR  Wings  Due  to  Circulation 
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Figure  1.2:  Nonlinear  Lift  of  LAR  Wings  Due  to  Wingtip  and  Leading  Edge  Vortices 


1.3  Scope  of  Present  Work 

1.3.1  Aerodynamic  Data  and  Analysis 

The  first  section  of  this  report  provides  a  summary  of  wind  tunnel  data  corresponding  to 
wings  of  low  aspect  ratio  and  various  planforms  at  low  Reynolds  numbers.  This  set  of  data 
includes  aspect  ratios  between  0.5  and  2.0,  four  distinct  wing  planforms,  and  Reynolds  numbers 
in  the  range  of  70,000  to  200,000.  The  results  consists  of  lift,  drag,  and  pitching  moment 
coefficients  as  a  function  of  angle  of  attack.  Analysis  of  this  experimental  data  reveals  important 
characteristics  of  LAR  wings  at  low  Re,  including  definite  nonlinearities,  high  stall  angles  of 
attack,  and  a  great  influence  of  aspect  ratio  and  wing  planform  on  aerodynamic  performance. 


1.3.2  Optimization  Code  and  Analysis  Procedure 

The  second  section  of  this  report  presents  one  of  many  possible  real-world  applications  of 
the  experiment  results  for  the  current  research.  The  data  of  force  and  moment  coefficients  versus 
angle  of  attack  for  each  of  the  planforms  and  aspect  ratios  is  implemented  into  an  aerodynamic 
analysis  program,  which  itself  is  part  of  a  design-optimization  procedure  applicable  to  small 
unmanned  aerial  vehicles  and  micro  air  vehicles.  The  definition  of  an  example  design  problem  is 
presented  first,  as  well  as  identification  of  design  variables  and  aircraft  performance  measures  of 
merit.  A  design-optimization  procedure  based  on  genetic  algorithms  is  discussed  next.  Finally, 
the  aerodynamic  analysis  code  that  determines  the  performance  of  a  given  set  of  design  variables 
is  presented.  This  analysis  procedure  makes  use  of  both  empirical  relations  based  on  wind  tunnel 
data  and  also  applies  conventional  aerodynamic  analysis  tools. 
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Chapter  2 

Experimental  Setup 


2.1  Wing  Models 

The  main  objective  of  the  present  study  is  to  determine  the  aerodynamic  characteristics  of 
LAR  wings  as  influenced  by  the  following  three  principal  factors: 

1.  Wing  planform. 

2.  Aspect  ratio. 

3.  Reynolds  number. 

In  order  to  study  the  effect  of  wing  planform,  four  distinct  geometries  were  chosen.  The  planforms 
used  were  rectangular,  Zimmerman,  inverse  Zimmerman,  and  elliptical.  The  Zimmerman  and 
inverse  Zimmerman  planforms  are  based  on  the  wingtip  shapes  designed  by  Zimmerman  in  the 
early  1930  s  (Zimmerman  [15])  and  consist  of  two  half-ellipses  joined  at  either  the  quarter- 
root-chord  location  (in  the  case  of  the  Zimmerman  planform)  or  at  the  three-quarter-root-chord 
location  (in  the  case  of  the  inverse  Zimmerman).  Figure  2.1  shows  an  illustration  of  how  the 
Zimmerman  planform  geometries  are  generated.  In  this  diagram,  two  half-ellipses  are  joined  at 
the  quarter-root-chord  location.  One  ellipse  has  semi-major  axis  a3  and  semi-minor  axis  ai  while 
the  other  has  semi-major  axis  a2  and  semi-minor  axis  a3. 


Figure  2.1:  Formation  of  the  Zimmerman  Planforms 
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For  each  of  the  four  planforms,  a  total  of  seven  aspect  ratio  values  were  tested  varying  from 
0.50  to  2.00  in  increments  of  0.25.  Figure  2.2  shows  the  planform  of  the  wings  tested  in  wind 
tunnel  experiments. 


Figure  2.2:  Wing  Planforms  Used  in  Experiments 


All  models  have  a  thickness-to-chord  ratio  of  1.96%  and  5-to-l  elliptical  leading  and  trailing 
edges  (see  Figure  2.3).  The  rectangular  planforms  have  flat  side  edges  and  elliptical  leading  and 
ti  ailing  edges  while  the  non-rectangular  planforms  have  elliptical  edges  around  the  entire  model. 

Foi  the  rectangular  wing  of  aspect  ratio  2.00,  two  additional  leading  edge  shapes  were  studied. 
These  were  a  10°  and  a  20°  tapered  leading  edge  shape,  as  shown  in  Figure  2.4  and  2.5.  Each 
of  these  tapered  shapes  is  actually  two  distinct  L.E.  geometries  depending  on  the  orientation  of 
the  wing  in  the  tunnel.  More  details  of  how  orientation  yields  different  L.E.  flow  geometries  are 
given  in  Section  3.6.  The  trailing  edges  of  these  modified  models  are  still  5-to-l  elliptical  in 


Figure  2.3:  Leading  Edge  Geometry:  5-to-l  Elliptical  L.E. 
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Figure  2.4:  Leading  Edge  Geometry:  10°  Tapered 
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Figure  2.5:  Leading  Edge  Geometry:  20°  Tapered 


Finally,  two  fuselage  body  models  were  used  to  determine  the  increase  in  drag  coefficient 
geneiated  by  the  piesence  of  a  fuselage.  The  fuselages  were  designed  to  be  similar  in  shape 
and  fineness  ratio  to  the  type  of  fuselages  usually  found  in  MAVs  and  small  UAVs.  The  two 
models  have  fineness  ratios  of  7.0  and  8.5.  Both  models  have  a  square  cross-section  with  rounded 
corners,  and  a  nose  cone  which  simulates  the  shape  of  a  typical  MAV  front  section.  Figure  2.6 
shows  a  picture  of  the  fuselage  bodies  attached  to  rectangular  and  inverse  Zimmerman  wings. 


Figure  2.6:  Fuselage  Bodies 
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2.2  Wind  Tbnnel 


Wind  tunnel  experiments  were  performed  using  a  low-turbulence,  indraft,  open-circuit  wind 
tunnel  located  at  the  University  of  Notre  Dame’s  Hessert  Center  for  Aerospace  Research.  This 
wind  tunnel  has  a  contraction  ratio  of  20.6  to  1,  a  rectangular  inlet  contraction  cone  designed  for 
low  turbulence  intensities,  and  a  2  ft  by  2  ft  test  section  (61  cm  by  61  cm).  Over  the  speed  range 
used  for  most  experiments,  the  turbulence  intensity  in  the  test  section  has  been  determined  to  be 
less  than  0.05%.  More  details  of  this  wind  tunnel  are  available  in  Mueller  [16].  A  schematic  of 
the  tunnel  is  shown  in  Figure  2.7. 


Figure  2.7:  Schematic  of  the  Low-Speed  Wind  Tunnel 

2.3  Force  Balance  and  Sting  Balance 

2.3.1  Platform  Force  Balance 

All  lift  and  diag  measurements  were  made  using  a  specially  designed  two-component  platform 
force  balance.  The  balance  is  mounted  externally  on  top  of  the  wind  tunnel  test  section.  A  sting 
and  sting  arm  connect  the  balance  to  the  model.  The  forces  are  measured  through  two  independent 
platforms  (one  for  lift  and  one  for  drag).  The  platforms  are  attached  to  two  flexures  on  which 
strain  gages  are  mounted  in  a  full  Wheatstone  bridge  configuration.  The  flexures  are  sized  to 
allow  a  resolution  of  1  gram,  which  allows  the  balance  to  measure  lift  and  drag  forces  of  wings 
models  down  to  a  Reynolds  number  of  approximately  20,000.  A  computer-controlled  stepper 
motor  and  associated  gearing  system  is  used  for  changing  the  angle  of  attack  of  the  model. 
Figure  2.8  shows  a  picture  of  the  platform  force  balance. 

The  models  were  mounted  vertically  and  attached  to  the  force  balance  sting  through  a  sting 
arm  in  the  trailing  edge  of  the  model  (see  Figure  2.9).  A  streamlined  sting  covering  protected 
the  sting  from  the  freestream  velocity,  such  that  the  only  forces  experienced  by  the  force  balance 
were  the  forces  on  the  model. 
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Figure  2.8:  Platform  Force  Balance 


Figure  2.9:  Mounting  of  Wing  Models  with  Platform  Force  Balance 


2.3.2  Five-Component  Sting  Balance 

A  five-component  internal  sting  balance  made  by  Aerotech  ATE  was  used  to  measure  normal 
force  and  pitching  moment.  This  balance  has  strain  gages  arranged  in  full-Wheatstone-bridge 
configurations  and  is  capable  of  measuring  normal  and  side  forces  as  well  as  pitching,  rolling, 
and  yawing  moments.  A  picture  of  the  balance  and  its  mounting  sleeve  (used  to  secure  the  wings 
to  the  balance)  is  shown  in  Figure  2.10. 


Figure  2.10:  Aerotech  Five-Component  Sting  Balance  and  Mounting  Sleeve 

The  wing  models  were  mounted  vertically  and  attached  directly  to  the  sting  balance  mounting 
sleeve.  A  machined  nose  cone  was  attached  to  the  front  of  the  sleeve  to  minimize  the  disturbance 
to  the  flow.  A  computer-controlled  stepper  motor  was  used  to  change  the  angle  of  attack  of 
the  wings.  Figure  2.11  shows  a  typical  mounting  arrangement  for  measurements  with  the  sting 
balance  (Zimmerman  model  of  aspect  ratio  2.00  is  pictured). 


Figure  2.11:  Mounting  of  Wing  Models  with  Sting  Balance 
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2.4  Signal  Conditioning  and  Data  Acquisition 

The  voltage  output  from  each  of  the  strain  gages  of  either  the  platform  force  balance  or 
the  sting  balance  was  processed  using  a  strain-gage  conditioning,  low-noise  amplifier  with  a 
10  Hz  low-pass  filter.  Typical  gains  used  for  the  platform  force  balance  were  in  the  order  of 
2,000  while  for  the  sting  balance,  the  gains  were  closer  to  4,000.  All  data  was  acquired  using  a 
PC-based  data  acquisition  system  running  the  LABVIEW®  5  graphical  programming  language. 
A  United  Electronics  Industries  12-bit  A/D  card  was  used.  Throughout  the  experiments,  8,000 
data  samples  were  acquired  quasi-simultaneously  per  input  channel  at  a  frequency  of  1000  Hz. 
Both  mean  and  standard  deviation  were  calculated  from  these  8,000  samples  per  channel.  The 
standard  deviation  is  used  in  the  uncertainty  analysis. 


2.5  Experimental  Procedures 

The  angle  of  attack  of  the  wing  models  was  varied  from  a  =  -10°  to  a  large  positive  angle 
between  a  =  30°  and  a  —  50°,  in  increments  of  1°.  Smaller  increments  of  angle  of  attack  were 
used  near  a  =  0°.  The  wing  was  then  brought  back  to  a  =  0°  to  determine  whether  hysteresis 
was  present  or  not.  No  hysteresis  was  observed  in  any  of  the  measurements.  At  each  angle  of 
attack,  the  voltage  of  an  electronic  manometer  connected  to  a  pitot-static  probe  in  the  tunnel  was 
also  measured  in  order  to  determine  the  dynamic  pressure  to  be  used  for  non-dimensionalization 
of  the  forces  and  moments. 

For  all  of  the  experiments  involving  the  force  balance  or  the  sting  balance,  an  offset  and  drift 
measurement  was  made  at  the  beginning  and  end  of  all  experiments.  The  drift  was  assumed  to 
vary  linearly  throughout  the  experiment. 

The  chord  Reynolds  numbers  presented  are  all  nominal  values  (not  corrected).  For  the 
experiments  of  lift,  drag,  and  moment  versus  angle  of  attack,  the  freestream  velocity  in  the 
test  section  was  always  adjusted  with  the  model  at  a  =  0°  to  yield  the  desired  nominal  Reynolds 
number.  The  force  and  moment  coefficients  presented  in  this  report  have  all  been  corrected  for 
wind  tunnel  blockage  (solid  blockage,  wake  blockage  and  streamline  curvature)  according  to  the 
techniques  presented  by  Pankhurst  and  Holder  [17].  Their  techniques  are  equivalent  to  techniques 
described  by  Rae  and  Pope  [18]. 
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2.6  Uncertainty 

The  two  main  sources  of  uncertainty  for  the  force  and  moment  coefficients  are  the  quantiza¬ 
tion  error  of  the  data  acquisition  card  and  the  uncertainty  arising  from  the  standard  deviation  of 
a  given  mean  strain  gage  output  voltage.  The  Kline-McClintock  [19]  technique  for  error  propa¬ 
gation  was  used  to  evaluate  all  uncertainties  in  the  aerodynamic  coefficients.  This  procedure  can 
yield  very  large  percentage  uncertainties  for  small  coefficients  approaching  zero.  For  larger  coef¬ 
ficients  (at  a  >  10°),  these  percentage  uncertainties  are  in  the  order  of  6%  for  CL,  CD  and  CM- 
As  an  example  of  typical  error  bounds  for  the  wind  tunnel  data  presented  in  this  report,  curves 
of  CL,  CD,  and  CM  (along  with  error  bars  corresponding  to  the  95%  uncertainty  bounds  of  each 
data  point)  are  plotted  in  Figures  2.12,  2.13,  and  2.14  versus  angle  of  attack.  The  rectangular 
wing  of  AR  =1.00  at  a  Reynolds  number  of  100,000  is  used  as  the  example  since  it  represents 
typical  values  for  the  worse-case  conditions  of  uncertainty.  As  the  Reynolds  number  or  aspect 
ratio  increases,  the  relative  uncertainty  values  decrease  due  to  the  larger  forces  involved.  The 
uncertainty  in  the  angle  of  attack  was  determined  to  be  of  the  order  of  0.5°  -  0.7°  for  most  test 
cases. 


Figure  2.12:  Uncertainty  in  CL,  Rectangular  Wing  of  AR  =1.00,  Re  =  100,000 
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Figure  2.13:  Uncertainty  in  CD,  Rectangular  Wing  of  AR  =1.00,  Re  =  100,000 


Figure  2.14:  Uncertainty  in  CM,  Rectangular  Wing  of  AR  =1.00,  Re  =  100, 000 
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Chapter  3 


Analysis  of  Experimental  Data 


3.1  Discussion  of  Aerodynamic  Results 

Appendix  A  contains  plots  of  the  experimental  results.  The  data  is  arranged  into  two  cate¬ 
gories:  grouped  by  aspect  ratio  and  grouped  by  wing  planform.  In  this  manner,  the  first  grouping 
is  appropriate  for  determining  the  effect  of  planform  for  wings  of  the  same  aspect  ratio  while  the 
second  grouping  is  appropriate  for  determining  the  effect  of  aspect  ratio  for  the  same  planform. 
Some  observations  about  the  data  are: 

1.  Wings  of  aspect  ratio  below  1.25  have  highly  nonlinear  lift  curves  that  are  characterized 
by  high  stall  angles  of  attack  and  non-constant  lift-curve  slopes  in  the  pre-stall  region. 

2.  Above  an  aspect  ratio  of  1.25,  most  of  the  planforms  exhibit  lift  curves  that  are  more  linear. 
As  predicted  by  theory,  the  higher  the  aspect  ratio,  the  more  linear  the  relationship  between 
lift  and  angle  of  attack. 

3.  For  aspect  ratios  below  1.50,  a  definite  performance  advantage  is  seen  for  the  inverse 
Zimmerman  and  rectangular  planforms  over  the  elliptical  and  Zimmerman  planforms.  This 
is  particularly  noticeable  when  comparing  the  Cl  versus  Cd  plots  for  constant  aspect  ratio 
(Figures  A.3,  A.7,  A. 13,  and  A. 19).  This  observation  is  further  discussed  in  Section  3.1.1. 

4.  Above  an  aspect  ratio  of  1.50,  the  difference  between  planforms  becomes  less  clear  and  is 
only  truly  apparent  at  very  high  angles  of  attack  beyond  stall. 

5.  Aspect  ratio  is  by  far  the  most  important  parameter  affecting  the  aerodynamic  characteristics 
of  LAR  wings  at  low  Re.  Wing  planform  is  the  next  most  important  factor,  followed  by 
Reynolds  number.  Leading  edge  shape  (discussed  in  more  detail  in  Section  3.6)  is  also 
critical.  However,  the  present  research  has  not  included  a  comprehensive  study  of  the  effect 
of  leading  edge  shape. 
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3.1.1  Lift  Performance  Related  to  Tip  Vortices 

The  effect  of  wing  planform  on  the  lift  characteristics  of  LAR  wings  for  a  given  aspect 
ratio  and  Reynolds  number  can  be  further  analyzed  by  introducing  the  parameter  xmax  span. 

max  span  is  the  chordwise  location  (measured  from  the  leading  edge)  of  maximum  wingspan, 
non-dimensionalized  by  the  root  chord  of  the  model.  For  the  Zimmerman,  elliptical,  and  inverse 
Zimmerman  wings,  xmax  span  is  0.25,  0.50,  and  0.75,  respectively.  For  the  rectangular  wings, 
Zmax  span  is  taken  to  be  1.0  rather  than  0.  In  essence,  xmax  span  provides  an  indirect  measure  of 
the  distance  between  the  wingtip  vortices  as  they  develop  over  the  wing  and  travel  downstream. 
It  was  determined  from  flow  visualization  experiments  that  the  distance  between  the  wingtip 
vortices  varied  proportionally  with  the  chordwise  location  of  maximum  span. 

This  can  best  be  described  by  use  of  the  sketch  and  photographs  of  Figure  3.1.  For 
wing  shapes  in  which  the  maximum  span  is  located  upstream  of  the  half-root-chord  location 
O^max  span  5;  0.5),  the  tip  vortices  are  seen  to  first  develop  at  the  location  of  maximum  wingspan. 
The  vortices  then  follow  the  outline  of  the  wing  up  to  a  point  and  separate  from  the  wing.  In 
contrast,  for  wings  with  xmax  span  greater  than  0.5,  the  vortices  are  seen  to  separate  from  the 
wing  at  the  location  of  maximum  span.  Thus  the  vortices  of  wings  with  xraax  span  >  0  5  are 
further  apart  than  those  of  wings  with  :rmax  span  <  0.5.  Figure  3.2  shows  another  set  of  flow 
visualization  photographs  which  outline  the  shape  and  relative  distance  of  the  wingtip  vortices 
over  low  AR  wings. 

Hoemer  [13]  suggested  that  the  lift  performance  of  low  aspect  ratio  wings  improves  as  the 
distance  between  the  wingtip  vortices  increases.  One  direct  way  to  measure  lift  performance  is 
to  compare  the  lift-curve  slope,  CLa,  for  each  wing  and  aspect  ratio.  The  lift-curve  slope  is 
calculated  through  a  least-squares  linear  regression  procedure  applied  to  the  experimental  data 
for  -10°  <  a  <  10°.  For  this  limited  range  of  angles  of  attack,  the  lift  curve  is  “linear  enough” 
to  allow  a  valid  determination  of  CLa.  As  can  be  seen  from  the  representative  case  at  a  Reynolds 
number  of  100,000  (Figure  3.3),  CLa  is  seen  to  increase  as  xmax  span  increases.  A  more  detailed 
analysis  of  the  lift  and  drag  characteristics  of  LAR  wings  can  be  found  in  Torres  and  Mueller  [20]. 
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Figure  3.1:  Separation  of  Wingtip  Vortices  Schematic  and  Comparison  with  Flow  Visualization 
of  Three  Planforms  (AR  =  1.0,  Re  w  70, 000,  a  =  15°) 
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Figure  3.2:  Flow  Visualization  of  Wingtip  Vortices  (AR  =  1.0,  Re  «  70,000,  a  =  15°) 
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Figure  3.3:  Lift  Curve  Slope,  CLa  vs.  xmax  span  for  Re  =  100, 000 
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3.2  Maximum  Lift  Coefficient  and  Stall  Angle  of  Attack 

An  important  characteristic  of  LAR  wing  aerodynamics  is  their  high  stall  angles  of  attack  and 
relatively  high  maximum  lift  coefficients.  The  experimental  data  of  Appendix  A  is  used  in  this 
section  to  define  CLmax  and  astaii  for  all  the  planforms  and  all  aspect  ratios  studied  in  this  work. 
Before  discussing  the  results,  some  issues  regarding  the  selection  of  the  maximum  lift  coefficient 
m  experiment  results  should  be  addressed.  In  most  cases,  the  maximum  lift  coefficient  can  easily 
be  determined  as  the  point  when  a  clear  stall  of  the  wing  is  observed.  This  is  particularly  true 
for  the  wings  of  higher  aspect  ratios.  However,  the  lower  aspect  ratio  wings  are  characterized 
by  highly  nonlinear  lift  curves.  In  some  cases,  a  point  of  inflection  in  the  lift  curve  exists  and 
the  lift  coefficient  decreases  slightly  before  increasing  again,  as  exemplified  in  Figure  3  4  for 
the  rectangular  planform  of  AR  =0.5.  In  other  cases,  the  point  of  inflection  is  not  present,  and 
the  lift  coefficient  increases  continuously  until  stall  at  a  very  high  angle  of  attack  as  seen  in  the 
inverse  Zimmerman  case  of  AR  =0.5  shown  in  Figure  3.5. 
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Figure  3.4:  Lift  Curve  of  Rectangular  Wing,  AR  =0.5,  Re  =  100,000 
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Figure  3.5:  Lift  Curve  of  Inverse  Zimmerman  Wing,  AR  =0.5,  Re  =  100, 000 
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In  an  effort  to  be  consistent,  the  maximum  lift  coefficient  is  taken  to  be  the  point  at  which 
the  lift  coefficient  reaches  a  local  maximum  and  then  decreases  in  value,  even  if  it  increases  in 
value  at  higher  angles  of  attack.  If  no  decrease  in  CL  is  observed,  such  as  in  Figure  3.5,  the 
^mai  is  taken  to  be  the  maximum  value  of  Cl  (approximately  1.25  in  the  example  of  Fig  3  5) 
It  is  suggested  that  Figures  A.41,  A.47,  A.53,  and  A.59  be  reviewed  as  references  to  what  is 
considered  Cxmax  in  each  case. 

Figures  3.6  and  3.7  plot  the  values  of  Cxmax  and  astaii  for  each  of  the  four  planforms 
as  a  function  of  aspect  ratio.  It  is  evident  from  these  graphs  that  wing  planform  has  significant 
influence  in  the  maximum  lift  characteristics  of  LAR  wings.  Specially  noticeable  is  the  difference 
between  the  rectangular  and  ellipse-based  planforms. 


AR 

Figure  3.6:  CLmax  of  Four  Wing  Planforms  and  All  Aspect  Ratios  (Re  =  100, 000) 
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Figure  3.7:  Stall  Angle  of  Four  Wing  Planforms  and  All  Aspect  Ratios  (Re  —  100, 000) 
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3.3  Aerodynamic  Center 

3.3.1  Background 

One  key  aspect  of  low  aspect  ratio  wing  aerodynamics  is  the  location  of  the  aerodynamic 
center  as  a  function  of  wing  planform,  aspect  ratio,  leading  edge  geometry,  Reynolds  number, 
and  angle  of  attack.  The  AC  location  yields  useful  information  about  the  structure  of  the  flow 
around  the  wing,  specially  in  the  nonlinear  lift  regime. 

The  location  of  the  AC  can  be  calculated  using  sting  balance  data  of  normal  force  and 
pitching  moment  taken  about  the  c/ 4  location  of  each  wing.  Before  continuing,  however,  a  few 
definitions  should  be  stated: 

Center  of  Pressure:  The  center  of  pressure  (xCP  or  hCP  =  *f-)  is  defined  as  the  point  along 
the  chordwise  direction  of  the  wing  at  which  the  lift  and  drag  forces  can  be  assumed  to 
be  acting.  In  other  words,  the  center  of  pressure  is  the  location  along  the  wing  where 
the  pitching  moment  is  identically  equal  to  zero.  For  symmetrical  airfoils,  the  center  of 
pressure  is  expected  to  lie  close  to  the  quarter-chord  location  (with  respect  to  the  wing’s  c) 
and  furthermore  it  is  expected  to  be  independent  of  angle  of  attack.  For  a  cambered  airfoil 
the  center  of  pressure  can  shift  significantly  as  the  angle  of  attack  changes. 

Aerodynamic  Center:  The  aerodynamic  center  (xAC  or  hAC  =  ^f1)  is  defined  as  the  point 
along  the  chordwise  direction  of  the  wing  at  which  the  pitching  moment  is  independent  of 
angle  of  attack,  that  is,  =0. 

Neutral  Point:  The  neutral  point  is  similar  to  the  aerodynamic  center  except  that  it  applies  to 
an  entire  aircraft.  That  is,  the  neutral  point  is  the  location  along  the  length  of  an  aircraft 
at  which  the  pitching  moment  is  independent  of  angle  of  attack. 

Note  that  for  a  wing  with  no  twist  and  with  airfoil  sections  which  are  symmetric,  the  center 
of  pressure  and  the  aerodynamic  center  are  one  and  the  same.  The  assumption  of  an  aerodynamic 
center  that  is  independent  of  angle  of  attack  and  which  lies  near  c/4  is  generally  found  to  be 

valid  for  wings  of  high  aspect  ratio.  As  the  aspect  ratio  decreases,  hAc  can  no  longer  be  assumed 
to  be  independent  of  a. 


3.3.2  Methods  for  Calculating  the  Location  of  the  Aerodynamic  Center 

Results  of  wind  tunnel  experiments  using  the  5-component  sting  balance  can  be  used  to 
determine  the  location  of  the  aerodynamic  center.  The  applicable  experimental  data  consists  of 
normal  force  and  pitching  moment  measured  at  the  quarter-chord  location  with  respect  to  the 
mean  aerodynamic  chord  of  the  wing.  Figure  3.8  shows  a  schematic  of  the  forces  and  moments 
involved  in  AC  analysis. 
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Figure  3.8:  Aerodynamic  Center  Analysis 

Transferring  the  forces  and  moments  from  the  quarter-chord  to  a  point  x^c  located  upstream 
of  the  quarter-chord  location  yields  the  following  expression  for  the  pitching  moment  about  the 
AC: 


Mac  =  Me/ 4  -  (zg/4  -  xAC)  N.  (3.1) 

Dividing  both  sides  of  this  equation  by  \pV^Sc  leads  to  the  nondimensional  form  of 
Eqn.  3.1: 


CWac  —  Cm  —  (0.25  —  hAc)  Cm,  (3.2) 

where  hAc  =  and  where  the  subscript  c/4  has  been  dropped  from  the  moment  coefficient 
for  simplicity.  There  are  several  ways  in  which  the  value  of  hAC  can  determined.  The  most 
straightforward  method  uses  the  knowledge  that  for  a  wing  with  no  camber  or  twist,  the  aero¬ 
dynamic  center  and  the  center  of  pressure  (or  center  of  lift)  are  one  and  the  same.  Furthermore, 
the  moment  at  the  AC  is  identically  equal  to  zero.  Setting  Eqn.  3.2  equal  to  zero  and  solving  for 
hAQ  yields 


hA  c  =  0.25-^.  (3.3) 

This  derivation  is  by  far  the  simplest  but  it  has  the  disadvantage  that  for  angles  of  attack  close 
to  zero,  the  value  of  Cm  goes  to  zero.  Therefore,  use  of  this  equation  yields  a  1/x  relationship 
as  a  goes  to  zero.  It  is  expected  to  be  accurate  only  for  moderate  to  high  angles  of  attack  where 
the  values  of  Cm  and  Cm  are  relatively  large. 
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Another  method  for  determining  the  location  of  the  aerodynamic  center  is  to  differentiate 
Eqn.  3.2  with  respect  to  angle  of  attack  while  assuming  that  hAC  is  a  constant  with  respect  to  a. 
By  definition  of  the  aerodynamic  center,  -Cq*c  =  0.  Thus, 


da 


QCm 

da 


-  (0.25  -  hAC)  ^ 


=  0. 


Solving  Eqn.  3.4  for  hAc  yields 


(3.4) 


hAc  =  0.25  - 


Cm a 
CNa 


(3.5) 


This  equation  does  not  have  a  discontinuity  near  a  =  0°  because  the  slope  of  the  normal  coefficient 
curve  is  not  zero  at  zero  angle  of  attack.  If  the  assumption  that  hAc  is  not  a  function  of  a  is 
correct,  this  equation  should  yield  a  constant  value  for  hAC  at  all  angles  of  attack  before  stall. 
Figure  3.9  shows  plots  of  hAC  versus  a  obtained  using  Eqns.  3.3  and  3.5  for  the  rectangular  wing 
of  AR  =1.0  at  Re  =  140,000.  The  values  of  C^a  and  Cm0  were  obtained  by  representing  the 
curves  of  Cm  and  Cm  versus  a  as  polynomials  of  order  4  and  10  respectively.  The  slope  was 
then  calculated  by  differentiating  the  polynomials  once  with  respect  to  a.  More  details  of  this 
polynomial  slope  method  are  presented  in  Section  3.3.3.  As  can  be  seen  in  the  graph  of  Figure  3.9, 
neither  equation  for  hAc  is  accurate  throughout  the  whole  range  of  angles  of  attack.  The  curve 
obtained  using  Eqn.  3.3  has  a  large  discontinuity  for  a  <  3°,  while  the  curve  corresponding  to 
Eqn.  3.5  is  definitely  not  constant  within  the  region  of  angles  of  attack  before  stall.  Therefore, 
it  can  be  concluded  that  the  assumption  of  hAc  being  independent  of  a  is  not  accurate  for  LAR 
wings  at  low  Re. 


Figure  3.9:  Calculation  of  hA c  using  Eqns.  3.3  and  3.5,  AR  =  1.0  Rectangular  Planform 
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Yet  another  approach  that  takes  into  account  h ac’s  dependence  on  a  is  to  differentiate 
Eqn.  3.2  with  respect  to  a  while  assuming  that  hAQ  is  not  a  constant.  This  yields 


®Cmac 

da 


—  -(0.25  -ftAC)^+C„a,*A0 


da 


da 


da 


=  0. 


(3.6) 


This  equation  can  be  solved  for  ^4^: 


= sb [(0'26  - hkc)  Cn-  -  -  e.7> 

At  a  —  0  ,  G\r  =  0  and  the  term  in  Eqn.  3.6  goes  to  zero.  Thus  a  forward  integration 
scheme  can  be  used  that  uses  the  initial  condition  at  a  =  0°: 


[^Ac] 


c*=0c 


=  hfl  =  0.25 


l^Ma  _  &={ 


a=0° 


[CkUoo 


0.25  - 


and  the  following  first-order  forward  Euler  method  to  integrate  the  angle  of  attack: 


(3.8) 
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*  =  0,1, 2, 3, 


(3.9) 


Figure  3.10  plots  the  values  of  hAC  obtained  through  the  use  of  Eqn.  3.9  as  well  as  Eqn.  3.3. 
It  can  be  concluded  from  this  graph  that  the  derivation  of  Eqn.  3.9  is  accurate  for  low  angles  of 
attack  and  the  one  of  Eqn.  3.3  is  accurate  at  higher  angles  of  attack.  There  is  a  region  in  which 
both  methods  overlap. 


Figure  3.10.  Calculation  of  hA c  using  Eqns.  3.3  and  3.9,  AR  —  1.0  Rectangular  Planform 
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3.3.3  Final  Method  Used  for  Calculating  hAc 

The  discussion  of  the  previous  section  provides  some  insight  into  the  complications  associated 
with  the  calculation  of  the  aerodynamic  center  for  LAR  wings  at  low  Re.  The  method  used  for 
this  report  to  determine  hAC  is  a  combination  of  Eqns.  3.3  and  3.9.  The  procedure  is  as  follows: 

1.  For  a  given  wing  planform,  aspect  ratio,  and  Reynolds  number,  wind  tunnel  data  exists  for 
CN  and  CM  as  a  function  of  angle  of  attack.  The  range  of  a  is  typically  between  -10° 
and  40°  with  increments  of  0.5°  for  -5°  <  a  <  10°  and  of  1°  for  all  other  angles. 

2.  A  polynomial  of  order  4  is  fitted  in  a  least  squares  sense  to  the  data  of  C A  vs.  a  for 
-3°  <a  <  15°.  Similarly,  a  polynomial  of  order  12  is  fitted  to  the  data  of  CM  vs.  a 
for  -3  <  a  <  15°.  The  reason  for  using  polynomials  instead  of  actual  experimental 
data  is  that  the  natural  scatter  of  the  data  generates  changes  in  slope  which  are  artificial 
in  nature.  The  polynomial  approximations  of  the  data  provide  a  smoothed  curve  and  ease 
the  task  of  calculating  slopes  since  all  that  is  required  is  to  differentiate  the  polynomial 
once  with  respect  to  a.  As  an  example,  Figure  3.11  shows  plots  of  CN  and  CM  and 
their  corresponding  polynomial  approximations  for  the  rectangular  planform  of  AR  =1  0 
at  Re  =  140, 000. 


Figure  3.11:  Polynomial  Approximations  to  CN  and  CM,  AR  =  1.0  Rectangular  Planform 


3.  The  value  of  hAC  at  a  =  0°  is  calculated  using  Eqn.  3.8.  The  values  of  the  slopes  of  CN 
and  Cm  are  obtained  by  differentiation  of  the  polynomials  obtained  in  the  previous  step. 
The  slopes  are  evaluated  at  a  =  0°. 

4.  Eqn.  3.9  is  used  to  integrate  forward  in  a.  The  integration  is  carried  out  for  0°  <  a  <  15° 
in  increments  of  0.1°. 

5.  For  angles  of  attack  greater  than  15°,  Eqn.  3.3  is  used.  The  values  of  CN  and  CM  are 
taken  directly  from  the  actual  experimental  data  (not  the  polynomial  approximation). 
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6.  For  most  of  the  wings  tested,  the  transition  from  the  differential  approach  of  Eqn.  3  9  to 
the  direct  approach  of  Eqn.  3.3  at  a  =  15°  is  seamless.  In  some  cases,  discontinuities  do 
occur,  but  they  are  in  general  small  and  do  not  significantly  affect  the  overall  trend  of  the 
curve  of  hAQ  versus  a. 

Figure  3.12  below  plots  the  final  result  of  this  procedure  for  the  rectangular  wing  of  AR  =1.0 
at  Re  =  140, 000.  The  transition  at  a  =  15°  is  relatively  smooth  and  the  value  of  hAC  is 
considered  valid  for  all  angles  of  attack  up  to  stall.  All  the  graphs  of  hAC  versus  a  shown  at  the 
end  of  this  section  were  obtained  using  the  method  outlined  above. 


Figure  3.12:  Calculation  of  hAC  Using  a  Combination  of  Eqns.  3.3  and  3.9 

The  uncertainty  in  the  value  of  hAc  is  difficult  to  calculate  analytically  as  the  procedure 
applied  to  determine  the  location  of  the  AC  involves  polynomial  curve  fitting,  forward  Euler 
integration,  and  combinations  of  more  than  one  technique.  In  lieu  of  an  analytical  uncertainty 
analysis,  a  computational  one  was  used.  That  is,  the  computer  program  which  implements  the 
procedure  outlined  in  this  section  was  run  repeatedly  (1000  times).  For  each  iteration,  the  values 
of  the  CN,  CM  and  a  were  randomly  varied  by  their  known  uncertainties.  Thus  1000  graphs  of 
hAc  vs.  a  were  obtained.  The  mean  and  standard  deviation  of  these  curves  provides  an  estimate 
of  the  uncertainty  of  the  value  of  hAc.  The  rectangular  wing  of  AR  =1.0  at  Re  =  100,000  was 
used  as  a  typical  test  case.  The  uncertainty  in  hAC  was  found  to  vary  between  2%  for  angles  of 
attack  close  to  0°  and  0.2%  for  angles  of  attack  above  15°.  In  general,  the  maximum  uncertainty 
of  hAC  in  the  graphs  presented  in  the  next  few  pages  is  no  larger  than  approximately  0.005  near 
a  =  °°  and  can  be  as  small  as  0.0005  for  a  >  15°.  Figures  3.13  through  3.22  at  the  end  of 
this  section  plot  the  value  of  hAC  versus  angle  of  attack  for  all  four  planforms  and  six  aspect 
ratios.  The  graphs  group  the  data  into  either  constant  aspect  ratio  (useful  to  determine  the  effect 
of  planform)  or  constant  planform  (useful  to  determine  the  effect  of  aspect  ratio). 
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3.3.4  Discussion  of  Aerodynamic  Center  Results 

The  plots  on  the  following  pages  provide  great  insight  into  the  lift-generating  mechanisms 
which  operate  in  the  regime  of  LAR  wings  at  low  Reynolds  number.  The  most  striking  observation 
that  should  be  made  is  the  drastic  shift  in  the  location  of  the  AC  as  angle  of  attack  increases. 
As  an  example,  consider  the  rectangular  wing  of  AR  =1.00  shown  in  Figure  3.18.  For  this 
wing,  hAC  is  close  to  0.17  at  a  =  0°  but  increases  to  0.40  at  a  =  40°.  Since  for  uncambered 
wings,  the  aerodynamic  center  and  the  center  of  lift  are  one  and  the  same,  it  can  be  concluded 
that  the  chordwise  location  at  which  the  lift  force  is  acting  shifts  downstream  as  the  angle  of 
attack  increases.  This  is  an  expected  trend  based  on  the  linear-nonlinear  theory  of  LAR  wings 
(see  Section  1.2  for  details).  For  low  angles  of  attack,  most  of  the  lift  generated  by  the  wing  is 
in  the  form  of  circulation  lift,  as  seen  in  wings  of  high  aspect  ratio.  Therefore,  it  is  expected 
that  the  location  of  the  AC  will  be  close  to  25%  mark  at  these  angles.  However,  as  the  angle 
of  attack  increases,  the  lift  mechanism  shifts  to  that  generated  by  the  wingtip  vortices.  These 
vortices  generate  lift  by  creating  low-pressure  sections  on  the  upper  surface  of  the  wing,  usually 
more  on  the  downstream  end  of  the  wing  (the  wingtip  vortices  increase  in  size  and  strength  as 
they  travel  downstream  along  the  upper  surface  of  the  wing).  It  would  be  expected  that  because 
the  lift  generated  at  higher  angles  of  attack  is  due  to  wingtip  vortices  (which  are  stronger  near  the 
trailing  edge),  the  location  of  the  center  of  lift  would  shift  downstream.  The  amount  by  which 
the  AC  shifts  towards  the  trailing  edge  is  a  direct  indication  of  how  much  nonlinear  lift  is  being 
generated  by  the  wing.  Wings  of  AR  =0.75  have  a  much  steeper  shift  of  AC  with  a  than  wings 
of  AR  =2.00.  If  one  were  to  plot  /iAC  versus  a  for  a  wing  of  AR  =10.00,  this  curve  would  be 
expected  to  be  a  horizontal  line  centered  close  to  0.25.  hA c  would  not  be  expected  to  change 
until  astaii  is  reached. 

Another  important  conclusion  that  can  be  drawn  from  the  AC  analysis  relates  to  the  effect 
of  planform  and  how  some  planforms  are  more  “linear”  than  others.  Zimmerman  and  elliptical 
planforms  are  seen  to  have  a  value  of  hAC  which  is  generally  lower  than  that  of  rectangular  and 
inverse  Zimmerman  planforms.  This  is  probably  due  to  the  fact  that,  as  discussed  in  Section  3.1.1, 
zim  and  ell  planforms  have  a  weaker  wingtip  vortex  system  than  the  one  found  in  ziminv  and 
red  planforms.  As  such,  it  would  be  expected  that  the  AC  of  the  planforms  with  the  weaker 
vortices  would  be  farther  forward  than  the  AC  of  wings  with  strong  vortices.  This  is  indeed  the 
case. 
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Figure  3.16:  hA0  vs.  a,  AR  =1.50,  Rec  =140,000 
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Figure  3.19:  /iac  vs.  a,  rectangular  planform,  Rec  =140,000 


a  (degrees) 

Figure  3.20:  /iac  vs.  a,  Zimmerman  planform,  Rec  =140,000 
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Figure  3.21:  h^c  vs.  a,  Inverse  Zimmerman  planform,  Rec  =140,000 
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Figure  3.22:  /iac  vs.  a,  Elliptical  planform,  Rec  =140,000 
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3.4  Effect  of  Reynolds  Number 

The  effect  of  varying  Reynolds  number  on  the  lift  and  drag  coefficients  of  LAR  wings  at 
low  Re  are  generally  small.  In  most  cases,  the  most  significant  effect  of  Reynolds  number  is 
associated  with  the  maximum  lift  coefficient  and  stall  angle  of  attack.  Figures  3.23  to  3.26 
plot  results  corresponding  to  the  rectangular  planform  of  AR  0.5,  1.0,  and  2.0  at  two  different 
values  of  Re.  The  trends  observed  for  these  particular  cases  are  similar  to  those  seen  in  all  other 
planforms  and  aspect  ratios.  It  can  be  seen  that  both  lift  and  drag  coefficients  increase  as  the 
Reynolds  number  increases.  Analysis  of  Figure  3.25  reveals  that  except  for  angles  of  attack  near 
stall,  the  effect  of  Reynolds  number  is  within  uncertainty  bounds  of  the  data.  This  figure  also 
shows  most  clearly  the  increase  in  stall  angle  of  attack,  specially  on  the  lower  aspect  ratio  cases. 


Figure  3.23:  Cl  vs.  a,  Rectangular  Planform,  Effect  of  Re 


3.5  Effects  of  Camber 

A  series  of  experiments  were  also  conducted  that  studied  the  effect  of  camber  on  LAR  wings. 
Brown  [21]  provides  a  more  detailed  set  of  results  than  available  in  this  report.  Figures  3.27 
through  3.30  reveal  some  general  trends.  In  these  experiments,  some  of  the  planform  shapes 
presented  in  Chapter  2  were  modified  to  have  4%  circular  camber.  Based  on  the  results  shown 
in  the  figures  of  this  section,  the  following  conclusions  about  the  effect  of  camber  can  be  stated: 
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Figure  3.24:  Cd  vs.  a.  Rectangular  Planform,  Effect  of  Re 

1.  Wings  with  camber  had  higher  lift  and  lower  drag,  yielding  significantly  more  efficient 
lifting  configurations  than  those  of  their  uncambered  counterparts 

2.  The  angle  of  attack  at  stall  and  the  maximum  lift  coefficient  were  seen  to  increase  dramat¬ 
ically  in  the  case  of  cambered  wings 

3.  The  most  significant  improvement  of  aerodynamic  efficiency  by  means  of  addition  of  cam¬ 
ber  was  seen  to  occur  for  the  wings  of  aspect  ratio  2.00.  The  performance  improvement 
(in  terms  of  lift-drag  polar)  of  wings  of  lower  AR  was  not  as  large. 

In  general,  all  experimental  results  point  to  the  conclusion  that  LAR  wings  at  low  Re  benefit 
from  camber  just  as  much  as  wings  of  high  aspect  ratio  at  higher  operating  Reynolds  numbers. 
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3.6  Effects  of  Leading  Edge  Shape 

A  preliminary  study  of  the  effects  of  leading  edge  shape  was  conducted.  The  objective  of 
this  study  was  to  determine  whether  leading  edge  shape  should  be  looked  at  more  closely  in 
future  research  projects.  This  analysis  also  reveals  how  sensitive  a  MAV’s  performance  is  to  the 
shape  of  the  leading  edge  and  how  much  care  must  be  taken  in  the  manufacture  of  wings.  The 
experimental  results  presented  in  this  section  correspond  to  a  rectangular  wing  of  AR  =2.00  at 
Reynolds  numbers  of  100,000  and  200,000.  The  wing  has  either  of  three  leading  edge  shapes: 
5-to-l  elliptical,  10°  wedge,  or  20°  wedge.  Depending  on  the  orientation  of  the  wing  (right  side 
up  or  upside  down),  each  of  the  wedge  leading  edge  shapes  is  actually  two  distinct  geometries 
(see  Figure  3.31). 
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20°  wedge  “down” 


Figure  3.31:  Leading  Edge  Geometries 

Figures  3.32  and  3.33  shows  plots  of  the  normal  force  coefficient  and  pitching  moment 
coefficient  about  the  quarter-chord  as  a  function  of  angle  of  attack  for  the  AR  2.0  rectangular 
planform.  The  effect  on  CV  is  not  very  large.  In  many  cases,  differences  in  normal  force 
coefficient  between  one  leading  edge  and  another  lie  within  uncertainty  bounds  of  the  experimental 
data.  The  most  significant  change  in  CN  is  visible  at  stall.  The  stall  angle  and  maximum  normal 
force  coefficient  for  the  wings  with  “down”  leading  edges  is  much  lower  than  that  of  the  wings 
with  “up”  leading  edges.  This  is  expected  as  the  “down”  leading  edges  form  an  extremely  sharp 
geometry  on  the  upper  surface  of  the  wing.  As  such,  flow  separation  occurs  earlier. 

In  contrast,  the  effect  on  the  moment  coefficient  is  quite  large  at  all  angles  of  attack.  Partic¬ 
ularly  interesting  in  Figure  3.33  is  the  shift  of  the  angle  of  attack  at  which  the  moment  coefficient 
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begins  to  decrease.  We  will  define  an  angle  of  attack  ctpeak  as  the  a  at  which  the  moment  coef¬ 
ficient  reaches  a  maximum.  As  an  example,  for  the  10°  “down”  leading  edge  at  Re  =  100,000, 
«peak  =  4.5°.  For  both  “down”  leading  edges,  CM  has  an  apeak  near  a  =  4.5°.  For  the  10° 
up  leading  edge,  apeak  is  near  a  —  8.5°  while  for  the  20°  “up”  leading  edge,  apeak  is  closer  to 
®  ^  *s  ProPosed  that  the  reason  for  this  shift  in  apeak  is  due  to  flow  separation.  For  the 

“up”  leading  edges,  an  adverse  pressure  gradient  (in  the  geometric  sense)  on  the  upper  surface 
of  the  wing  does  not  exist  until  the  angle  of  attack  goes  beyond  the  angle  of  the  leading  edge 
wedge.  Thus,  it  is  expected  that  the  flow  will  remain  attached  to  the  upper  surface  up  until  that 
angle.  For  the  “down”  leading  edges,  separation  is  expected  to  occur  immediately  because  the 
wedge  angle  is  essentially  zero.  Evidence  for  the  proposed  explanation  can  be  seen  in  the  curves 
of  Figure  3.33.  The  correlation  is  not  perfect,  as  exemplified  by  the  fact  that  apeak  does  not  occur 
at  20°  for  the  20°  wedge.  Nevertheless,  a  trend  definitely  exists  that  links  the  angle  of  inclination 
of  the  leading  edge’s  upper  surface  and  the  value  of  c*peak. 

Another  interesting  result  of  Figure  3.33  is  the  effect  of  Reynolds  number  on  apeak.  For  all 
the  geometries  shown  in  this  figure  (except  for  the  10°  “down”  one),  a  jump  in  the  CM  vs.  a 
curve  occurs  as  the  Reynolds  number  increases  from  100,000  to  200,000  (for  example,  this  jump 
occurs  at  a  =  1.75°  for  the  elliptical  leading  edge). 

The  observations  and  proposed  explanations  provided  in  this  section  are  preliminary.  Not 
enough  experimental  data  is  available  to  throughly  explain  the  effects  of  leading  edge  shape.  One 
sure  conclusion  of  this  phase  of  the  project  is  that  further  study  of  leading  edge  influence  on 
aerodynamic  characteristics  should  be  made  in  the  future.  Such  a  study  is  beyond  the  scope  of 
the  present  work.  Future  research  at  the  Hessert  Center  for  Aerospace  Research  could  focus  on 
this  aspect  of  LAR  wing  aerodynamics. 
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a  (degrees) 

Figure  3.32:  Effect  of  Leading  Edge  Shape  on  C^,  Rectangular  Planform  AR  =2.0 


Chapter  4 


Design-Optimization  Part  I:  Design  Variables  and 

Optimization  Algorithm 


4.1  Introduction  to  MAV  Design-Optimization 

The  previous  chapters  of  this  report  have  summarized  wind  tunnel  data  for  wings  with  the 
following  criteria: 

1.  Four  distinct  planforms  (red,  zim,  ziminv,  ell). 

2.  Seven  values  of  aspect  ratio  (0.50,  0.75,  1.00,  1.25,  1.50,  1.75,  and  2.00). 

3.  Three  values  of  Reynolds  numbers  (70,000,  100,000,  and  140,000). 

4.  Two  values  of  camber  (0%  and  4%). 

5.  Three  fuselage-body  configurations  (no  fuselage-body  and  two  rounded  bodies  with  differ¬ 
ent  fineness  ratios). 

The  data  consists  mostly  of  curves  of  Cl,  Cm,  Cd,  and  Cm  versus  angle  of  attack  for  each  wing 
planform,  aspect  ratio,  and  Reynolds  number  combination.  The  location  of  the  aerodynamic 

center  of  any  of  the  wings  tested  as  a  function  of  angle  of  attack  can  be  determined  using  this 
same  data.  c 

It  is  important  to  stress  that  this  wind  tunnel  data  is  directly  applicable  to  real  world  aero¬ 
dynamic  design  problems.  The  data  itself  does  not  yield  much  practical  information  unless  it 
is  incorporated  into  some  type  of  useful  aerodynamic  analysis  tool.  This  is  the  objective  of  the 
next  two  chapters.  In  this  chapter,  an  example  design  challenge  typical  to  MAV  applications  is 
introduced.  After  appropriate  restrictions  are  imposed,  the  design  variables  and  aircraft  perfor¬ 
mance  characteristics  are  defined.  The  design  problem  is  well  suited  to  the  use  of  optimization 
algorithms.  Therefore,  a  genetic  algorithm  global  optimization  technique  is  suggested  to  solve 
the  problem.  The  next  chapter  shows  how  the  experimental  data  obtained  in  this  work  can  be 
incorporated  into  an  aerodynamic  analysis  procedure  applicable  to  MAVs.  This  analysis  tool  is 
used  in  the  optimization  program  in  order  to  generate  optimum  designs. 

The  application  presented  in  the  next  two  chapters  is  by  no  means  meant  to  be  a  final  product 
or  to  represent  the  only  possible  use  of  the  wind  tunnel  data.  It  is  meant  to  show  how  useful 
engineering  decision-making  tools  can  be  extracted  from  wind  tunnel  data. 
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4.2  Design  Problem  Statement 

4.2.1  Summary  of  Design  Optimization  Problem 

The  design  problem  considered  in  this  chapter  incorporates  requirements  and  mission  objec¬ 
tives  that  are  common  to  many  current  MAV  implementations.  Even  though  mission  objectives 
vary  greatly  depending  on  the  application  and  (sometimes  more  importantly)  on  the  end  user, 
many  MAV  missions  share  common  characteristics.  In  the  discussions  that  follow,  the  design 
problem  is  defined  to  be  the  following: 

The  designei  is  presented  with  the  challenge  of  developing  the  smallest  remote-controlled 
air  vehicle  that: 

a)  can  carry  a  wireless  video  transmitter/video  camera  payload, 

b)  has  an  endurance  of  at  least  10  minutes  and  a  surveillance  range  of  at  least  0.8 
Km  (0.5  mi)  and 

c)  is  maneuverable  enough  to  be  able  to  capture  an  image  of  a  target  on  the  ground. 

The  size  of  the  vehicle  is  defined  explicitly  as  the  maximum  linear  dimension  between  any 
two  points  located  on  the  aircraft  while  it  is  airborne. 

4.2.2  Detailed  Mission  Requirements 

The  mission  requirements  listed  above  can  be  further  analyzed  in  order  to  obtain  more 
quantitative  design  objectives.  These  are: 

Video  Payload:  The  MAV  must  be  able  to  cany  video  components  that  have  a  range  of  at  least 
0.8  Km  (0.5  mi).  Video  components  that  satisfy  these  specifications  can  be  selected  from 
off-the-shelf  components  and  thus  their  dimensions  and  weights  are  known. 

Range/Endurance:  For  the  purposes  of  this  example  problem,  the  propulsion  system  of  the 
MAV  is  restricted  to  be  an  internal  combustion  engine.  The  particular  engine  was  chosen 
based  on  the  expected  takeoff  weight  of  the  MAV.  The  fuel  consumption  of  the  engine  is 
known  from  engine  tests;  therefore,  the  required  fuel  weight  needed  for  an  endurance  of 
10  minutes  can  be  calculated  directly. 

Stability  and  Control:  The  airplane  should  be  able  to  cruise  in  a  trimmed  flight  condition  (that 
is,  the  pitching  moment  about  the  center  of  gravity  should  be  close  to  zero).  Also,  the 
pitching  moment  derivative  with  respect  to  angle  of  attack  at  cruise  must  be  negative 
in  sign  (longitudinal  stability  requirement).  For  control,  the  airplane  must  carry  control 
surface  actuators  (servos),  a  remote-control  radio  receiver,  and  batteries  (which  power  both 
the  control  system  and  the  video  system).  These  components  are  also  predefined  and  their 
weights  and  dimensions  are  known. 
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4.2.3  Restrictions 


Certain  restrictions  must  be  placed  on  the  aircraft  in  order  to  apply  the  semi -empirical  analysis 

procedure  and  in  order  to  simplify  the  design  process.  These  restrictions  are: 

Fixed  Wing:  The  MAV  is  restricted  to  a  fixed-wing  aircraft  that  may  or  may  not  have  a  horizontal 
stabilizer.  The  results  section  of  this  chapter  discusses  the  complications  arising  in  tail-less 
airplane  configurations.  However,  for  the  present  work,  the  optimizer  is  allowed  to  select 
a  design  set  that  has  zero  horizontal  tail  area. 

Wing  and  Tail  Thickness,  Camber,  and  Edge  Shape:  The  main  aerodynamic  analysis  is  based 
on  the  interpolation  of  wind  tunnel  data  of  low  aspect  ratio,  flat-plate  wings  (details  of  this 
procedure  are  given  in  Section  5.1.1).  As  such,  the  wing  and  tail  surfaces  must  be  restricted 
to  zero-camber,  a  thickness  to  chord  ratio  of  1.96%,  and  5-to-l  elliptical  leading  and  trailing 
edges  so  that  they  match  the  parameters  of  the  wing  models  used  in  the  experiments. 

Vertical  Location  of  Tail:  It  is  assumed  that  if  a  tail  is  present,  the  tail  is  located  behind  the 
wing  (not  a  canard  configuration).  In  addition,  the  tail  is  restricted  to  be  located  below  the 
wing.  The  vertical  separation  between  the  wing  and  tail  is  set  equal  to  the  fuselage  height. 
The  reason  for  this  restriction  is  that  flight-testing  of  a  variety  of  MAV  configurations  has 
revealed  that  significantly  more  stable  flight  can  be  achieved  if  the  horizontal  tail  is  located 
well  clear  of  the  wing  wake. 

Tractor  Configuration:  The  engine  is  assumed  to  be  at  the  nose  of  the  airplane  in  a  tractor 
configuration. 

Predefined  Components:  The  main  components  of  the  MAV  (engine,  propeller,  fuel  tank,  bat¬ 
teries,  video  transmitter,  video  camera,  RC  receiver,  and  RC  servos)  have  been  selected 
based  on  commercially  available  products  that  are  as  small  and  light  as  possible  but  still 
satisfy  mission  requirements.  This  restriction  is  particularly  important  in  terms  of  the  engine 
because  in  traditional  airplane  design  procedures,  the  selection  of  the  propulsion  system  is 
part  of  the  design  process. 

Fixed  Take-off  Weight:  Since  the  components  to  be  carried  by  the  MAV  are  predefined,  the  size 
of  the  fuselage  that  will  house  the  components  is  also  essentially  fixed.  The  fuselage  is  the 
largest  contributor  of  structural  weight  of  an  MAV;  therefore  it  is  possible  to  approximately 
define  a  fixed  take-off  weight  that  includes  the  component  weights  as  well  as  estimates 
for  the  structural  weights  of  the  fuselage,  wing,  tail,  vertical  stabilizer,  and  miscellaneous 
accessories  (linkages,  hinges,  switches,  wiring,  etc.).  The  fixed  value  of  takeoff  weight  W 
for  the  MAV  to  be  designed  in  this  report  is  1.0  Newton  (0.22  lb),  equivalent  to  a  mass  of 
approximately  100  g  (3.52  oz). 
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4.3  Design  Variables,  Design  Parameters,  and  Aircraft  States 

4.3.1  Design  Variables 

For  an  MAV  with  the  configuration  limitations  mentioned  above,  the  design  variables  are 
defined  as  follows: 


Variable 

Symbol 

Type 

Main  Wing 

Wing  Planform 

Pw 

Discrete 

Aspect  Ratio 

ARw 

Continuous 

Planform  Area 

Sw 

Continuous 

Incidence  Angle 

iw 

Continuous 

Horizontal  Tail 

Tail  Planform 

PT 

Discrete 

Aspect  Ratio 

ARrp 

Continuous 

Planform  Area 

St 

Continuous 

Incidence  Angle 

%t 

Continuous 

Position  Variables 

Location  of  Tail 

xt 

Continuous 

Location  of  Propeller 

prop 

Continuous 

Aircraft  Angle  of  Attacl 

c 

Aircraft  Angle  of  Attack 

a 

Continuous 

Table  4.1:  Design  Variables 


Wing  and  tail  variables  which  are  treated  as  constants  in  this  example  are  the  thickness  to 
chord  ratio  (1.96%),  percent  camber  (0%),  and  leading  and  trailing  edge  shape  (5-to-l  elliptical 
shape).  As  mentioned  before,  the  reason  for  fixing  these  four  variables  as  constant  is  that  wind 
tunnel  data  is  used  in  the  primary  aerodynamic  analysis  through  an  interpolation  procedure.  In 
order  to  match  the  experiment  data  with  the  design  variables  of  the  MAV,  these  four  parameters 
are  set  equal  to  those  of  the  experiments. 

Note  that  the  planform  of  the  wing  and  tail  is  defined  as  being  a  discrete  variable.  This  is 
due  to  the  fact  that  the  wind  tunnel  data  used  in  the  main  aerodynamic  analysis  includes  only  four 
wing  planforms.  The  shape  parameters  associated  with  these  planforms  are  not  evenly  spaced 
and  thus  do  not  offer  a  simple  way  of  fitting  a  response  surface  based  on  wing  shape.  Therefore 
the  wing  and  tail  are  allowed  only  four  possible  shapes.  Nevertheless,  the  wing  area  and  aspect 
ratio  are  truly  continuous  variables  and  their  aerodynamic  performance  can  be  interpolated. 

The  location  variables  xt  and  ajpr0p  determine  the  location  of  the  horizontal  tail  and  propeller 
in  the  streamwise  direction.  These  two  variables  alone  define  the  placement  of  all  other  compo¬ 
nents  due  to  restrictions  placed  on  the  relative  location  of  components  inside  the  aircraft.  Section 
5.3  provides  details  of  these  restrictions  and  how  they  affect  the  aircraft’s  center  of  gravity. 

Finally,  the  vertical  position  of  the  wing  and  tail  is  not  defined  as  a  design  variable  but  is 
considered  a  parameter  instead.  See  Section  5.4.1  for  an  explanation  of  how  the  values  of  vertical 
distance  between  the  wing/tail  and  the  fuselage  centerline  ( zw  and  z?,  respectively)  are  assigned. 
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4.3.2  Design  Parameters 

Design  parameters  are  those  design  elements  that  cannot  be  directly  controlled  by  the  designer. 
For  the  example  design  problem  considered  in  this  report,  the  design  parameters  can  be  sub¬ 
classified  into  either  component  parameters  or  operational  parameters: 

Component  Parameters:  The  component  parameters  include  the  size,  weight,  and  center  of 
mass  of  the  components  to  be  carried  on  board  (engine,  propeller,  fuel  tank,  batteries, 
two  servos,  radio  control  receiver,  video  transmitter,  and  video  camera).  The  geometry 
and  mass-distribution  of  these  components  (given  in  Appendix  B)  are  inputs  to  the  overall 
design  process.  Furthermore,  the  total  takeoff  weight  of  the  aircraft  is  assumed  to  be 
constant  and  known  [equal  to  1.0  Newtons  (0.22  lb)]. 

Also  included  in  this  category  is  the  performance  of  the  engine  to  be  used  for  the  MAV. 
Experiments  with  this  engine  have  yielded  curves  of  thrust  versus  freestream  velocity.  This 
information  is  used  in  the  determination  of  the  cruise  speed. 

Finally,  the  weight  density  of  the  material  used  to  construct  the  wing  and  tail  is  also 
considered  a  component  parameter  and  is  used  in  center  of  gravity  calculations.  The  MAV 
prototypes  constructed  at  Notre  Dame  over  the  last  few  years  have  been  made  using  0.3175 
cm  (1/8  in)  foam  and  it  is  assumed  that  this  is  the  material  of  choice  for  the  design 
considered  in  this  report. 


Operational  Parameters:  The  principal  operational  parameters  are  air  density  and  freestream 
velocity,  which  are  needed  for  converting  non-dimensional  coefficients  into  forces.  Density 
is  a  user-supplied  value  while  freestream  velocity  is  calculated  iteratively  within  the  analysis 
program  as  described  in  Section  5.2.2.  Other  operational  parameters  could  be  air  viscosity 
or  wind  speed;  however  these  secondary  operational  parameters  are  not  considered  in  this 
report. 


4.3.3  Aircraft  States  (Performance) 

The  main  states  of  an  MAV  design  are  the  maximum  dimension,  the  total  lift  force  generated, 
and  the  overall  pitching  moment  and  trim  characteristics.  In  more  detail,  these  states  are: 

Maximum  Dimension:  The  maximum  dimension  of  a  given  design  depends  on  the  size,  aspect 
ratio,  and  shape  of  the  wing  and  tail,  the  location  of  these  two  surfaces,  and  the  location 
of  the  propeller  (from  which  the  position  of  all  other  components  is  derived).  Section  5.5 
discusses  the  details  of  the  determination  of  maximum  dimension. 

Total  Lift:  The  total  lift  generated  by  the  aircraft  is  dependent  on  the  size,  aspect  ratio,  shape, 
and  incidence  angle  of  the  wing  and  tail,  angle  of  attack  of  the  aircraft,  and  on  the  cruise 
velocity.  The  total  lift  force  acting  on  the  airplane  must  be  equal  to  the  takeoff  weight  of 
the  MAV,  W  (within  some  threshold).  In  addition,  the  angle  of  attack  of  the  wing  and 
tail  should  be  lower  than  the  stall  angle  of  attack  for  each  lifting  surface.  In  general,  the 
cruise  angle  of  attack  is  desired  to  be  much  smaller  than  the  stall  angle  of  attack  to  allow 
for  maneuvering.  This  additional  angle  is  defined  here  as  Asafety. 
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Pitching  Moment:  The  moment  about  the  center  of  gravity  must  be  equal  to  zero  (or  within  e 2 
of  zero)  for  steady  cruise  flight.  Also  related  to  pitching  moment  is  the  requirement  that 
the  derivative  of  pitching  moment  with  respect  to  angle  of  attack  be  negative  in  sign. 


4.4  Optimization  Problem 

4.4.1  Statement  of  Optimization  Problem 

The  discussion  in  the  previous  sections  provides  a  preliminary  basis  for  stating  the  complete 
design/optimization  problem  related  to  the  design  of  a  micro  air  vehicle  in  standard  optimization 
form.  The  constraints  have  been  converted  into  less  than  or  equal  constraints  in  order  to  facilitate 
the  use  of  optimization  algorithms. 

Objective  Function: 

Maximum  dimension  of  the  MAV,  Dmax  =$■  to  be  minimized. 

Design  Variables: 

x=[P\v>  ARw,  Sw>  iw,  Pt,  ARt,  St,  ir,  %t,  ^propjQl] 

Subject  to  the  Constraints: 

||  W  —  LtotalU  <  Cl 

OL  i\y  T  dfsafety  ®  stall  w  —  ® 

a  T  %T  T  (^safety  Q stall?  —  0 

ll^call  <  e2 

9Mcg  <  q 

da  ~ 

4.4.2  Design  Variable  Constraints 

The  design  variables  are  constrained  to  satisfy  the  following  requirements: 


Pwmin 

< 

Pw 

< 

PWmax 

ARWmin 

< 

ARw 

< 

AR'Wmax 

SWmin 

< 

Sw 

< 

SWmax 

iWmi  n 

< 

iw 

< 

i-Wmax 

PTmin 

< 

Pt 

< 

Pi  \nax 

ARWmin 

< 

ARt 

< 

ARTmax 

ST  . 

J- min 

< 

St 

< 

ST 

-L  max 

%T  • 

A  mm 

< 

It 

< 

^Tmax 

XTmir, 
-1  mm 

< 

xT 

< 

XTmax 

Xpropmin 

< 

xprop 

< 

Xpvopmax 

Qmin 

< 

a 

< 

&max 
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4.4.3  Genetic  Algorithm 

4.4.3.1  Summary  of  Genetic  Algorithm 

The  MAV  design-optimization  problem  described  thus  far  is  a  highly  nonlinear  design  prob¬ 
lem  with  many  design  variables.  Therefore  it  lends  itself  to  nonlinear  global  optimization  methods 
such  as  simulated  annealing  or  genetic  algorithms.  A  genetic  algorithm  (GA)  method  will  be  used 
in  this  work;  it  is  implemented  following  the  description  of  Belegundu  [22].  Genetic  algorithms 
are  applicable  for  the  maximization  of  a  function  / (a  ■ )  (where  x  is  the  vector  containing  the  de¬ 
sign  variables  X\,  x2, ...,  x n)  subject  to  upper  and  lower  bounds  on  x  (that  is,  xmin  <  x  <  xmax ). 
A  brief  description  of  the  algorithm  follows: 

1.  Each  continuous  design  variable  x{  is  represented  as  a  binary  number  of  m2  bits.  The 
interval  of  said  design  variable  (xmax.  -  xmin.)  is  therefore  partitioned  into  (2m2  -  1)  sub¬ 
intervals.  A  special  case  exists  for  the  two  discrete  variables  Pw  and  Pj*  (the  planforms 
of  the  wing  and  tail).  These  variables  can  only  take  four  values:  1,  2,  3,  or  4,  each 
representing  one  of  the  four  planforms  tested  in  wind  tunnel  experiments.  A  convenient 
way  to  represent  these  discrete  design  variables  is  to  define  the  number  of  bits  for  Pw  and 
Pt  to  be  equal  to  mi=2.  This  automatically  allows  only  four  possible  conditions.  The 
number  of  bits  of  the  other  continuous  variables  (m2)  is  still  a  parameter  at  the  control  of 
the  user. 

2.  A  population  of  Z  members  is  created  by  attaching  the  binary  representation  of  each  x{ 
end  to  end.  As  an  example,  one  member  of  a  population  with  m2  =  8  would  have  8 N  bits 
and  would  look  as  follows: 

101 1010101001 1 1 100110101, ,  .111Q1Q1Q 

X1  x2  X3  xN 

The  number  of  members  in  a  population  usually  varies  between  10  and  200.  Note  that  the 
example  above  ignores  the  complication  of  the  discrete  variables  having  a  different  number 
of  bits  than  the  number  of  bits  of  the  continuous  variables.  The  two  discrete  variables 
use  mi =2  bits  while  the  continuous  variables  use  m2- 8  bits.  For  this  case,  the  binary 
representation  of  a  population  member  would  be: 

v  10  J-  0 1 1 0 10  ljO  1 001 111001 10101^  11  .1001010T..1110101Q 

Xl  ^2  X3  X4  xs  Xq  XN 

Note  that  the  number  of  bits  in  a  member  of  the  population  for  the  discrete  variable  case 
is  Ambits  =  [2mi  +  m2(N  -  2)] 

3.  An  initial  population  is  selected  using  a  random  number  generator  that  randomly  selects  a 
1  or  a  0  for  each  of  the  iVbits  bits  in  each  member  of  the  population. 

4.  Each  member  of  the  population  is  converted  from  a  string  of  l’s  and  0’s  into  N  design 
variables  xi,  x2, ...,  xN  by  transforming  the  binary  representation  into  actual  values. 

5.  The  function  value  /  is  evaluated  for  each  set  of  design  variables  x  corresponding  to  each 
member  of  the  population. 
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6.  A  mating  pool  is  created  such  that  the  members  of  the  population  whose  function  value  is 
greater  have  a  greater  probability  of  “reproducing”.  This  is  done  as  follows: 

(a)  The  maximum  and  minimum  values  of  /  in  the  population  (/ft  and  fL)  are  used  to 
define  two  parameters  C  and  D: 

C  =  0.1/ft  -  1.1  f i  (4.1) 

D  =  max(l,fh  +  C)  (4.2) 

(b)  The  function  values  /*  for  each  member  of  the  population  are  then  scaled  using  C 
and  D : 

f'i  =  (4.3) 

(c)  The  parameter  S  is  defined  as  S  =  /•  for  i  =  1...Z. 

(d)  A  random  number  r  between  0  and  1  is  generated  Z  times.  Each  time,  the  index  j  is 
found  such  that 

E  fl<rS  <•£}',  (4.4) 

i=l  i=  1 

The  jth  member  of  the  population  is  copied  into  the  mating  pool. 

(e)  If  desired,  the  procedure  for  selecting  the  mating  pool  can  be  modified  such  that  the 
first  Abest  members  of  a  given  population  are  guaranteed  to  be  in  the  mating  pool. 
The  procedure  outlined  above  would  then  be  used  for  the  remaining  members  of  the 
mating  pool. 

7.  Crossover  is  performed  as  follows: 

(a)  Two  members  of  the  mating  pool  population  are  selected  to  be  “parents”. 

(b)  A  random  number  k  between  1  and  Nbits  is  generated. 

(c)  The  first  k  bits  of  the  two  parents  are  exchanged  and  they  become  the  first  k  bits 
of  the  “children”.  The  remaining  bits  of  the  children  are  copied  from  those  of  the 
parents.  As  an  example,  the  following  crossover  operation  would  take  place  if  k  is 
randomly  determined  to  be  13: 

parentl  100100101010101010001010... 

parents  P100100110111, 11110100101... 

- ^ - / 

child 1  010010011011101010001010... 
child2  /001001010101  11110100101... 

(d)  The  crossover  operation  is  performed  Z/2  times,  generating  a  new  child  population 
with  Z  members. 
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8.  Mutation  is  accomplished  by  pre-defining  a  mutation  probability  bp.  Every  time  that  a 
crossover  operation  is  performed,  a  random  number  rm  is  selected  between  0  and  1  for 
each  of  the  Ambits  bits  of  each  child.  If  rm  <  bp,  then  the  bit  is  changed  to  a  0  if  it  was  a  1 
or  to  a  1  if  it  was  a  0. 

9.  The  function  value  is  evaluated  for  each  member  of  the  new  child  population.  Steps  6 
through  8  are  repeated  AIqen  times  (equal  to  the  number  of  generations). 

4.4.3.2  Penalty  Functions 

The  MAV  optimization  problem  must  be  converted  from  a  constrained  minimization  problem 
to  a  form  that  is  applicable  to  the  GA  method.  This  requires  the  use  of  penalty  functions  to 
represent  the  constraints  of  the  problem.  The  five  constraints  (lift  constraint,  wing/tail  stall 
angle  constraints,  trimmed  flight  constraint,  and  longitudinal  stability  constraint)  are  converted 
to  penalty  functions  as  follows: 

1.  Lift  Constraint:  \\W  -  Ltota i||  <  ex 

A  penalty  function  P1  is  defined  such  that: 

PX  =  (W  -  Ltotal)2  (4.5) 

2.  Stall  Angle  of  Attack  Constraint:  a  +  ix  +  asafety  -  astaUx  <  0  (X  =  W  or  T) 

Penalty  functions  P2  and  P3  are  defined  such  that: 

P2  TTUlx(0 ,  (ct  T  %\\r  -\-  asafety  ^ stall  ) )  (4.6) 

P3  =  max( 0,  (a  +  iT  +  asafety  -  astaiiT))2  (4.7) 

3.  Trimmed  Flight  Constraint:  ||Mcg||  <  e2 
A  penalty  function  P4  is  defined  such  that: 

P4  =  (Mcg)2 *  (4.8) 

4.  Longitudinal  Stability  Constraint:  <  0 

A  penalty  function  P5  is  defined  such  that: 

P5  =  max(0,  (MCg&+A&  -  MCg& ))2  (4.9) 

where  MCGa+Aa  is  the  moment  about  the  CG  at  an  angle  of  attack  that  is  higher  than  the 
design  angle  by  an  (arbitrary)  amount  Act. 
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4.4.3.3  Statement  of  Optimization  Problem  in  GA  Form 

The  penalty  functions  described  in  the  previous  section  are  used  to  convert  the  unconstrained 
minimization  problem  of  Section  4.4.1  into  a  constrained  maximization  problem  which  is  in  the 
standard  genetic  algorithm  form: 

GA  objective  function: 

/($)  =  ~(D  max  +  T\P\  +  V2P2  +  T3P3  +  +  r5P5),  to  be  maximized. 

Subject  to: 

—min  —  S.  —  —max 

The  function  f(x)  is  calculated  for  a  given  set  of  design  variables  x  =  x1,x2,  ...,xN  using 
the  aerodynamic  analysis  procedure  described  in  Chapter  5.  The  factors  T\  through  r$  are  penalty 
function  multipliers  chosen  by  the  user.  The  values  of  the  penalty  function  multipliers  used  in  the 
optimization  for  this  report  are  n=15,000,  r2= 10,  r3=10,  74=100,000,  and  r5=l,000,000.  The 
reason  for  such  large  values  of  r4  and  r5  is  that  the  pitching  moment  about  the  CG  usually  has 
values  smaller  than  0.01  N  -  m.  The  corresponding  penalty  functions  must  be  scaled  such  that 
they  become  significant  in  comparison  to  the  maximum  dimension  or  the  other  penalty  functions. 

The  minimum  and  maximum  values  that  were  used  for  the  MAV  optimization  are  given  in 
Table  4.2. 


Variable 

Minimum 

Maximum 

Units 

Pw 

rect,  zim,  ziminv,  ell 

AR\y 

0.5 

2.0 

Sw 

25  (3.88) 

400  (62.00) 

cm2  (in2) 

iw 

0 

15 

degrees 

PT 

rect,  zim,  ziminv,  ell 

art 

0.5 

2.0 

St 

0 

200  (31.00) 

cm 2  (in2) 

it 

-15 

0 

degrees 

xT 

0 

15  (5.91) 

cm  (in) 

^prop 

0 

15  (5.91) 

cm  (in) 

a 

-5 

25 

degrees 

Table  4.2:  Minimum  and  Maximum  Values  for  Design  Variables 
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4.5  Results  and  Discussion 

4.5.1  Test  of  Genetic  Algorithm  Code 

A  FORTRAN  code  of  the  genetic  algorithm  was  written  which  implements  the  steps  of  the 
GA  code  as  described  in  Section  4.4.3.  This  code  was  tested  by  solving  an  optimization  problem 
that  has  a  known  solution.  An  objective  function  Ftest  was  defined  as: 


-Ftest  —  ~[(wi[Pw  -  Pw0 ])2  +  (w2[ARw  -  ARWo})2  +  ( w3[Sw  -  £V0])2  +  ~  iw0])2 

+  (w5 [Pt  ~  Pro})2  +  (w6[ARt  -  ARTq])2  +  (w7[ST  -  5r0])2  +  (u>8[iT  -  zTo])2 

T  ( Wq\xt  -^Tol)  P  (^lot^prop  ^propo])  T  (fUnjci  —  <5o])]2  (4.10) 

where  the  variables  with  the  subscript  0  are  constants  and  the  factors  wi  through  wn  are  constant 
multipliers  which  scale  each  term  so  that  it  contributes  approximately  equally  to  the  objective 
function.  The  maximum  value  of  Ftest  is  achieved  when  each  variable  is  equal  to  its  corresponding 
0-subscript  constant  (that  is,  Pw  —  P\y0>  ARw  =  ARw0,  Sw  =  *SW0,  iw  —  *wb>  etc.).  The  values 
of  the  constants  and  multipliers  used  in  the  test  case  are  shown  on  Table  4.3. 


Variable 

Target  (A0) 

GA  Result 

Scaling  multiplier  ( wx ) 

Pw 

1.000 

1.000 

1.000 

AR\y 

1.350 

1.356 

2.000 

Sw 

155.0 

155.1 

0.005 

iw 

7.500 

6.884 

0.200 

Pr 

3.000 

3.000 

1.000 

ARx 

0.850 

0.825 

2.000 

ST 

95.00 

102.4 

0.005 

it 

-4.500 

-4.790 

0.200 

•  Xt 

7.500 

6.940 

0.100 

^prop 

9.500 

8.673 

0.150 

a 

3.500 

3.413 

0.100 

Table  4.3:  Results  of  GA  Test  Case 


The  genetic  algorithm  parameters  used  were  a  population  of  50  members,  200  generations, 
12-bit  discretization  of  the  design  variables  (except  for  the  two  discrete  planform  variables  Pw 
and  Pt  which  use  2-bit  discretization),  and  a  mutation  probability  of  0.05.  The  lower  and  upper 
bounds  of  the  design  variables  were  the  same  as  those  in  Table  4.2.  Figure  4.1  shows  the 
convergence  history  of  the  test  case.  The  negative  of  the  objective  function  converges  to  a  value 
close  to  the  exact  optimum  as  the  number  of  generations  increases.  The  final  values  of  the  design 
variables  for  the  case  plotted  in  Figure  4.1  are  shown  on  Table  4.3.  Based  on  the  results  of  the 
test  case  for  the  GA  code,  it  can  be  concluded  that  the  genetic  algorithm  is  working  correctly.  It 
is  now  possible  to  apply  the  algorithm  to  the  actual  design  optimization  problem. 
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1 


0.9 


Generation  Number 

Figure  4.1:  Convergence  History  of  GA  Test  Case 

4.5.2  Preliminary  Results:  Problems  With  Discrete  Variables 

The  first  attempts  at  running  the  actual  optimization  revealed  some  shortcomings  of  the 
procedure  as  outlined  thus  far.  Convergence  almost  never  occured  regardless  of  the  values 
chosen  for  population  size,  number  of  generations,  mutation  probability,  or  bit-discretization.  In 
almost  all  cases,  the  function  value  fluctuated  wildly  and  did  not  seem  to  converge.  The  cause 
for  this  behavior  was  traced  back  to  the  two  discrete  variables  of  wing  and  tail  planform  (Pw  and 
PT).  As  mentioned  before,  these  variables  are  only  allowed  to  take  four  possible  values  which 
represent  each  of  the  four  planforms  tested  in  wind  tunnel  experiments.  The  problem  arises  from 
the  fact  that  the  aerodynamic  characteristics  of  two  planforms  differ  greatly  even  if  all  other 
factors  (aspect  ratio,  angle  of  incidence)  remain  constant.  The  GA  would  drive  all  continuous 
design  variables  to  an  optimum  for  a  given  combination  of  wing  and  tail  planform.  However,  if 
the  planform  variable  was  changed,  the  resulting  aircraft  would  differ  so  much  from  the  previous 
“near-optimum”  design  that  convergence  would  be  nearly  impossible. 

The  solution  to  this  problem  was  to  remove  the  two  discrete  variables  (wing  planform  and  tail 
planform)  from  the  list  of  design  variables  which  are  controlled  by  the  GA  code.  That  is,  for  each 
optimization  run,  the  planform  of  the  wing  and  tail  are  defined  and  kept  constant.  This  leaves  9 
continuous  design  variables  which  define  completely  the  MAV  design.  One  possibility  for  dealing 
with  the  remaining  two  discrete  planform  variables  is  to  run  16  cases,  each  case  corresponding  to 
a  combination  of  wing  planform  and  tail  planform.  The  smallest  maximum  dimension  of  the  16 
optimized  designs  would  be  considered  to  be  the  global  optimum.  For  the  purpose  of  this  report, 
solving  only  one  of  the  16  cases  proves  the  applicability  of  the  GA  optimization  code  to  the 
design  of  MAVs.  Performing  the  remaining  15  cases  would  then  be  a  straightforward  procedure. 
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It  was  decided  to  define  the  wing  to  be  of  inverse  Zimmerman  planform  and  the  tail  to  be  of 
rectangular  planform.  These  wing  shapes  are  kept  constant  and  the  remaining  9  design  variables 
are  used  in  the  GA  code  to  arrive  at  an  optimum  for  this  configuration. 

4.5.3  Results  of  Modified  Optimization  Problem 

4.5.3. 1  Design  Variables  of  Optimum  Design  (Converged  Solution) 

With  the  modified  genetic  algorithm  procedure  it  was  found  that  convergence  was  easier  to 
achieve.  Approximately  eight  out  of  every  ten  times  that  the  code  was  run  resulted  in  convergence 
to  a  feasible  solution.  There  are  a  number  of  factors  that  could  cause  lack  of  convergence  in 
the  unsuccessful  cases.  First  and  foremost,  the  objective  function  itself  is  highly  nonlinear. 
Small  changes  in  one  variable  have  major  impacts  on  the  performance  of  the  design.  It  becomes 
difficult  for  the  GA  to  minimize  the  maximum  dimension  of  the  aircraft  while  also  keeping  the 
penalty  functions  as  small  as  possible.  In  many  of  the  failure  cases,  the  optimizer  converged  to  a 
solution  which  was  not  feasible  (penalty  functions  active)  but  did  not  shift  from  this  solution.  A 
higher  value  of  the  mutation  probability  was  found  to  fix  this  problem  but  it  made  the  optimizer 
simply  jump  randomly  from  one  design  condition  to  another  one.  Other  factors  that  were  seen 
to  influence  convergence  were  the  penalty  function  multipliers  and  the  number  of  members  in 
a  given  population.  Some  more  experimentation  with  these  parameters  might  prove  helpful  in 
improving  the  chances  of  a  feasible  convergence  of  the  code. 

Overall,  converged  solutions  were  significantly  more  common.  Figure  4.2  shows  one  such 
successful  case.  It  plots  the  aircraft’s  maximum  dimension  versus  generation  number.  The  GA 
optimization  parameters  used  for  this  run  were  a  population  of  40  members,  60  generations,  12- 
bit  discretization  for  the  continuous  variables,  and  a  mutation  probability  of  0.06.  Table  4.41ists 
the  values  of  the  optimized  design  variables  (corresponding  to  the  fittest  member  of  generation 
number  60).  Figure  4.3  shows  a  planform  view  of  the  optimized  design. 

The  optimized  design  of  Figure  4.3  is  similar  to  MAV  prototypes  built  at  Notre  Dame. 
These  aircraft  have  successfully  flown  with  100  grams  takeoff  weight.  Specially  important  is 
the  fact  that  the  optimized  MAV  has  a  horizontal  tail.  There  were  other  optimization  cases  that 
converged  on  a  design  that  did  not  have  a  horizontal  tail.  That  is,  in  these  cases,  the  optimizer 
identified  a  tail-less  airplane  as  a  more  efficient  configuration  in  terms  of  maximum  dimension. 
This  is  logical  since  the  tail  is  not  actually  necessary  to  satisfy  the  very  limited  trimmed  flight 
constraint.  It  may  seem  that  the  tail  only  adds  to  the  maximum  dimension  without  providing 
any  additional  benefit.  However,  the  aerodynamic  analysis  used  as  part  of  the  GA  code  does  not 
include  a  very  important  effect  associated  with  tail-less  airplanes  with  low  aspect  ratio  (LAR) 
wings:  LAR  wings  usually  operate  at  relatively  large  angles  of  attack  (10°  in  some  cases).  At 
these  high  angles,  the  flow  separates  near  the  trailing  edge  of  the  wing  without  inducing  stall. 
Control  surfaces  (which  are  usually  located  at  the  trailing  edge  of  the  wing  in  tail-less  airplane 
configurations)  can  become  ineffective  when  immersed  in  this  separated  air  region.  Ineffective 
control  surfaces  yield  highly  unstable  flight  even  at  cruise  conditions.  Extensive  flight-testing  of 
MAV  prototypes  has  revealed  that  a  tail-less  configuration  can  sustain  stable  flight  for  a  takeoff 
weight  which  is  approximately  25%  smaller  than  a  tailed  configuration  of  the  same  maximum 
dimension.  Unfortunately,  at  the  moment  it  is  difficult  to  quantitatively  define  the  disadvantages 
of  tail-less  flight  within  the  aerodynamic  analysis  portion  of  the  code.  The  understanding  of 
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exactly  how  the  wing’s  wake  affects  the  control  of  the  MAV  is  not  very  vast  yet.  Therefore,  for 
the  present  version  of  the  optimization  algorithm,  the  problems  associated  with  tail-less  airplane 
configurations  are  ignored  and  a  tail-less  configuration  is  not  considered  in  any  way  inferior  to 
one  that  does  have  a  tail. 

The  design  defined  in  Table  4.4  is  considered  to  be  a  feasible  design  even  though  the  moment 
about  the  center  of  gravity  is  not  exactly  zero.  A  moment  of  0.0003  N  -mis  equivalent  to  a 
3  gram  weight  located  1  cm  from  the  center  of  gravity.  This  small  moment  is  not  expected  to 
change  the  trim  characteristics  of  the  MAV  considerably.  As  a  reference,  it  should  be  compared 
to  the  change  in  pitching  moment  arising  from  the  loss  of  fuel  weight  in  the  fuel  tank  as  a  flight 
progresses.  A  typical  location  for  the  fuel  tank  is  2  cm  forward  of  the  CG  (of  course,  the  exact 
location  of  the  tank  is  different  for  each  design).  The  fuel  tank  carries  7  grams  of  fuel.  Thus,  at 
the  end  of  the  flight,  a  change  in  pitching  moment  which  is  more  than  3  times  larger  than  0.0003 
N  -  mis  present  due  simply  to  fuel  usage.  It  is  therefore  considered  justifiable  to  treat  0.0003 
N  —  m  as  essentially  zero  moment. 
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Span  of  wing 


Span  of  tail 
Maximum  dimension 
Cruise  speed 
Lift 

_ Mcg 

Mcg“+a“ 


1.412 

103.0  (15.97) 
-9.68 

6.270  (2.47) 
3.950  (1.56) 
2.53 

20.72  (8.16) 


12.06  (4.75) 
21.20  (8.35) 
11.78  (26.36) 
1.024 
0.0003 
-0.0018 


cm  (in) 
cm  (in) 
m/s  ( MPH ) 
N 

N  —  m 
N  —  m 


Table  4.4:  Optimized  MAY  Design 


4.6  Conclusions  and  Recommendations  of  MAY  Optimization 


In  this  chapter,  a  genetic  algorithm  optimization  method  was  applied  to  the  design-optimization 
of  a  micro  air  vehicle.  Typical  mission  requirements  were  used  to  define  specific  constraints  and 
restrictions.  The  genetic  algorithm  code  was  found  to  converge  on  feasible  designs  that  met  all 
constraints  and  satisfied  all  mission  requirements  and  restrictions. 

In  some  cases,  convergence  of  the  algorithm  did  not  occur.  This  is  attributed  largely  to 
the  highly  nonlinear  character  of  the  objective  function  coupled  with  the  restrictive  nature  of 
the  constraints.  It  was  also  found  that  the  aerodynamic  performance  analysis  used  within  the 
GA  code  is  not  capable  of  modeling  the  unstable  flight  characteristics  that  can  occur  with  some 
tail-less  configurations.  This  is  an  important  aspect  of  the  design  because  tail-less  airplanes 
can  feasibly  represent  optimum  configurations  in  respect  to  maximum  dimension.  Yet  flying 
characteristics  of  certain  configurations  of  tail-less  aircraft  are  not  good.  Although  a  detailed 
study  of  aerodynamic  interactions  associated  with  tail-less  and  tailed  aircraft  configurations  is 
beyond  the  scope  of  this  work,  it  is  still  suggested  that  some  type  of  qualitative  analysis  be 
incorporated  in  an  effort  to  model  the  complications  of  tail-less  configurations  more  accurately. 
Flow  visualization  experiments  with  different  configurations  are  planned  for  the  near  future  and 
will  provide  qualitative  insight  into  the  flow  structure  of  separated  flow  regions  in  LAR  wings. 
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Chapter  5 

Design-Optimization  Part  II:  Aerodynamic  Analysis 


This  chapter  discusses  in  detail  the  procedure  used  to  calculate  a  value  for  the  objective 
function  used  within  the  optimization  scheme  based  on  a  set  of  design  variables  x.  The  analysis 
consists  of  the  following  steps: 

1.  Determine  the  values  of  CL  and  CD  of  the  wing  and  tail  given  each  surface’s  angle  of 
incidence,  aspect  ratio,  and  wing  planform. 

2.  Check  for  stall  condition  for  the  wing  and  tail. 

3.  Iteratively  determine  the  cruise  velocity  of  the  aircraft  and  the  lift  and  drag  forces  of  the 
wing  and  tail. 

4.  Calculate  the  location  of  the  center  of  gravity. 

5.  Calculate  the  pitching  moment  about  the  CG. 

6.  Check  for  longitudinal  stability  (negative  slope  of  Mqg  vs.  a  curve). 

7.  Determine  the  maximum  dimension  of  the  MAV  design. 

8.  Calculate  the  values  of  the  penalty  functions  based  on  violations  of  constraints  and  penalty 
function  multipliers. 

9.  Combine  the  maximum  dimension  and  penalty  functions  to  calculate  the  value  of  the  ob¬ 
jective  function  f(x)  for  the  given  design  defined  by  x. 


5.1  Interpolation  of  Aerodynamic  Coefficients  and  Stall  Angle 
Analysis 

The  first  step  of  the  procedure  (shown  schematically  in  Figure  5.1)  is  to  calculate  the  lift  and 
drag  coefficients  of  the  wing  and  tail  based  on  the  values  of  7  design  variables:  the  planform, 
aspect  ratio,  and  incidence  angle  of  the  wing  and  tail  and  the  angle  of  attack  of  the  aircraft.  The 
stall  angle  is  also  determined  as  a  function  of  these  variables  and  a  check  is  made  for  both  the 
wing  and  the  tail  to  check  for  a  stall  condition.  Note  that  the  stall  check  for  the  tail  does  not 
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include  the  effect  of  downwash.  The  change  in  effective  angle  of  inclination  of  the  local  flow 
of  the  horizontal  tail  due  to  downwash  is  expected  to  be  small.  Future  revisions  of  the  analysis 
procedure  will  include  the  downwash  effects. 


Check  for  stall 

[C*-H\V+  ^safety] ^  ^stnll  W? 


Wing  DV’s  (Pw,  ARW,  iw) 
Aircraft  DV’s  (a) 

Tail  DV’s  (PT,  ARt,  iT) 


► 


Interpolate  to  get 
C^w^nd  CDW 


Figure  5.1:  Analysis  Procedure: 


5.1.1  Determination  of  Lift  and  Drag  Coefficients  for  Wing  and  Tail 

Research  on  low  aspect  ratio  wing  aerodynamics  at  low  Reynolds  numbers  has  revealed  that 
analytical  expressions  for  lift  and  drag  coefficients  are  seldom  accurate.  Even  though  the  slope 
of  the  linear  region  of  the  lift  curve  can  be  approximated  quite  well,  most  LAR  wings  operate  in 
the  nonlinear  region  of  lift.  Analytical  methods  for  nonlinear  lift  have  been  found  not  to  compare 
well  with  wind  tunnel  experiment  data.  Correlations  of  experiment  drag  curves  with  analytical 
predictions  have  been  determined  to  be  particularly  inaccurate. 

A  proposed  solution  to  this  problem  is  to  use  the  wind  tunnel  data  presented  in  this  report 
as  the  basis  of  an  interpolation  scheme.  The  wind  tunnel  data  corresponds  to  flat-plate  wings 
of  aspect  ratios  of  0.5,  0.75,  1.0,  1.25,  1.5,  1.75,  and  2.0  and  four  distinct  wing  planforms  (rec¬ 
tangular,  Zimmerman,  inverse  Zimmerman,  and  elliptical).  This  experimental  data  consists  of 
curves  of  lift  and  drag  coefficients  versus  angle  of  attack  (-10°  <  a  <  40°  in  increments  of  1°). 
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The  chord-Reynolds  numbers  at  which  these  sets  of  data  were  taken  are  70,000  and  100,000, 
with  some  wings  having  been  tested  at  140,000.  Most  MAV  designs  with  a  takeoff  weight  of 
approximately  1  Newton  operate  at  Reynolds  numbers  between  100,000  and  200,000.  Experi¬ 
mental  data  has  revealed  little  change  of  the  aerodynamic  characteristics  from  Re  =  100, 000 
to  Re  =  140, 000.  Thus  the  lift  and  drag  coefficients  at  Re  =  100, 000  are  considered  as  the 
correct  values.  All  wing  models  tested  had  a  thickness  to  chord  ratio  of  1.96%,  camber  of  0%, 
and  5-to-l  elliptical  leading  and  trailing  edges. 

The  lift  and  drag  coefficients  of  a  particular  wing  or  tail  design  are  considered  to  be  a  function 
of  three  variables  only:  wing  planform,  aspect  ratio,  and  angle  of  inclination  with  respect  to  the 
incoming  flow  (aefr).  Since  the  planforms  are  discrete  variables,  Cl  and  Cd  are,  for  a  given 
planform,  actually  functions  of  only  two  independent  variables.  Note  however  that  the  angle  of 
inclination  of  the  wing  or  tail  with  respect  to  the  incoming  flow  is  itself  a  function  of  two  design 
variables:  a  and  iw  (or  a  and  iT  for  the  tail).  That  is,  aeffw  =iw  + a  and  aeffr  =  ir  +  a.  A 
linear  interpolation  procedure  is  used  to  obtain  the  value  of  lift  and  drag  coefficients  for  a  given 
set  of  design  variables.  The  interpolation  algorithm  is  the  following: 

1.  The  planform  of  the  wing  or  tail  is  read  and  the  database  of  wind  tunnel  data  corresponding 
to  that  planform  is  loaded. 

2.  The  lift  and  drag  coefficients  are  interpolated  linearly  at  the  desired  aeff  for  all  aspect 
ratios. 

3.  The  interpolated  coefficients  from  step  2  are  then  used  to  interpolate  the  final  value  of  CL 
and  Cd  at  the  desired  aspect  ratio.  Figures  5.2  and  5.2  show  examples  of  the  interpolation 
procedure. 

Once  the  Cl  and  Cd  are  known,  they  are  converted  to  dimensional  forces: 


Lw 

=  CLwl-pVjSw 

(5.1) 

Dw 

—  Cdw~pVoo2Sw 

(5.2) 

Ltj * 

=  Cl^pVJSt 

(5.3) 

Dt 

=  CDtv\pVoo2St , 

(5.4) 

where  rj  accounts  for  the  change  in  dynamic  pressure  at  the  tail  with  respect  to  the  dynamic 
pressure  at  the  wing.  In  the  first  implementation  of  this  design-optimization  procedure,  rj  is  set 
to  1.0. 

Voo  is  not  known  a  priori  and  must  be  calculated  in  an  iterative  manner.  See  Section  5.2.2 
for  details  on  how  this  iterative  procedure  is  carried  out. 
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5.1.2 


Determination  of  Stall  Angle  for  Wing  and  Tail 


The  angle  of  attack  at  stall  is  also  needed  within  the  optimization  program  In  Section  3  9 
graphs  of  astall  vs.  A*  for  each  of  the  four  planforms  were  presen”'  discuss  d  An 

Tat  o7r^:rhr,ar  ,o  ,he  °r  r in  ,he  ~  £ 

Z,  7  f  7taU;  u  ner’  a  Value  0f  astaI>  could  be  calculated  given  the  wing  planform 

undated  t  Ts  t  °f  Writing  ***  rep0rt’ the  aerodynamic  anflysis  code  had  not  been 
pckted  to  include  such  an  interpolation  procedure.  Instead,  the  average  value  (at  each  value  of 

AR)  for  all  four  w,ng  planforms  was  plotted  versus  AR  as  shown  in  Figure  5  4  An  emptaUv 

etata  taTr  •  °',0Se'y  approxinKtes  the  average  a"g'e  cufve  is  also  plotted.  The 
equation  ror  this  empmcal  approximation  is  F 


^staii  ~  -10  atan(4(AR  -  1.25))  +  28c 


p.oj 


OeJTTT  ‘he  interpola«™  technique,  Equation  5.5  is  used  within  the  analysis  program  to 
determme  the  maxmtum  allowable  angle  of  attack  of  the  wing  and  tail.  Future  vtstonfof  the 

unique  to  eaaeh  plaZn  d  dete™"ati™  °f  “*•»  «»»*  an  in.etpola.ion  method  that  is 


Figure  5.4:  Average  Stall  Angle  of  Attack  for  All  Planforms 
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5.2  Drag  Analysis  and  Cruise  Speed 

The  total  drag  force  is  determined  here  through  an  iteration  procedure.  The  cruise  velocity 
is  needed  to  calculate  the  total  drag  force,  and  in  turn,  the  thrust  available  from  the  propulsion 
source  (engine)  is  dependent  on  cruise  speed.  Thus,  in  an  iterative  manner,  the  cruise  speed 
at  which  the  drag  and  thrust  are  equal  to  each  other  is  calculated.  The  procedure  used  in  this 
analysis  is  diagrammed  in  Figure  5.5. 


L=CLjMaI(0.5p  VJ  Sw)  +  T  sin(a) 


Figure  5.5.  Analysis  Procedure:  Cruise  Speed  and  Lift  and  Drag  Forces 

5.2.1  Total  Drag  Force 

Thus  far,  the  drag  coefficients  of  the  wing  and  tail  are  known  based  on  the  interpolation 
procedure  of  Section  5.1.1.  An  additional  drag  increment  exists  due  to  the  presence  of  the 
fuselage,  engine,  and  vertical  tail,  as  well  as  due  to  aerodynamic  interactions  between  these 
items  and  the  wing/tail.  Because  the  fuselage  is  expected  to  have  the  greatest  influence  on 
the  total  airplane  drag  increment,  wind  tunnel  experiments  were  conducted  that  measured  the 
difference  in  drag  coefficient  between  a  wing  alone  and  a  wing  outfitted  with  a  typical  MAV 
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fuselage.  In  such  a  manner,  an  empirical  method  for  estimating  the  increase  in  drag  due  to  the 
fuselage  can  be  formulated.  The  experiments  were  done  using  rectangular,  Zimmerman,  and 
inverse  Zimmerman  wings  of  AR  1.0,  1.5,  and  2.0.  Figures  5.6  through  5.8  show  results  for  the 
rectangular  planforms  of  aspect  ratio  1.0, 1.5,  and  2.0.  In  each  plot,  two  curves  are  graphed.  One 
of  them  is  the  difference  in  drag  coefficient  between  the  wing  outfitted  with  a  fuselage  and  the 
wing  alone.  That  is,  this  is  the  increment  in  drag  coefficient  due  to  the  presence  of  the  fuselage. 
The  other  curve  plots  the  drag  coefficient  of  the  fuselage  by  itself  (without  a  wing).  The  drag 
coefficient  for  the  fuselage  is  nondimensionalized  by  the  same  reference  area  as  that  of  the  wing 
corresponding  to  each  graph.  It  can  be  seen  from  these  graphs  that  the  actual  increment  in  drag 
coefficient  due  to  the  fuselage  is  generally  greater  than  the  drag  of  the  fuselage  alone.  This  is 
likely  due  to  interference  effects  that  cannot  be  modeled  in  a  test  of  the  fuselage  alone. 


a  (degrees) 

Figure  5.6:  Increment  of  Drag  due  to  Fuselage,  Rectangular  Wing  of  AR  =1.0 


a  (degrees) 

Figure  5.7:  Increment  of  Drag  due  to  Fuselage,  Rectangular  Wing  of  AR  =1.5 
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Figure  5.8:  Increment  of  Drag  due  to  Fuselage,  Rectangular  Wing  of  AR  =2.0 


It  is  interesting  to  note,  though,  that  the  drag  increment  is  actually  lower  than  the  Cd  of  the 
fuselage  alone  for  the  wings  of  AR  =1.0  and  AR  =1.5  at  angles  of  attack  beyond  15°.  This 
is  again  expected  to  be  due  to  interference  effects  and  highly  nonlinear  aerodynamic  interactions 
occuring  between  the  wing  and  the  fuselage  body.  A  detailed  investigation  of  the  linear  and 
nonlinear  effects  associated  with  fuselage-body  drag  is  beyond  the  scope  of  the  current  research. 
Rather,  it  is  desired  to  obtain  an  approximation  to  the  additional  drag  generated  by  the  presence 
of  a  fuselage  body  on  a  lifting  wing.  Thus,  for  the  purpose  of  the  aerodynamic  analysis  presented 
herein,  the  average  increase  in  drag  coefficient  due  to  the  fuselage  is  taken  to  be  0.015  (based 
on  the  wing  reference  area).  This  additional  drag  coefficient  is  added  to  the  aircraft  total  drag 
budget,  yielding  a  total  drag  coefficient  for  the  entire  airplane. 


5.2.2  Thrust  Force  and  Cruise  Speed 

Once  the  total  aircraft  drag  coefficient  is  known,  the  drag  force  can  be  calculated  based  on 
the  freestream  velocity,  T4.  However,  14,  is  not  known  a  priori.  It  is  determined  through  an 
iterative  procedure  which  matches  the  thrust  from  the  engine  with  the  drag  force  of  the  airplane 
at  a  given  freestream  velocity.  Figure  5.9  shows  the  variation  of  engine  thrust  with  freestream 
velocity  (obtained  experimentally)  for  a  Cox  0.010  Tee  Dee  engine.  This  particular  engine  is 
typically  used  on  small  remote-controlled  airplanes  and  represents  the  most  adequate  off-the- 
shelf  propulsion  system  for  non-battery-powered  MAV  designs.  Again,  the  emphasis  here  is  on 
the  design-analysis  procedure,  not  on  the  details  of  each  step.  The  code  could  easily  be  changed 
such  as  to  accommodate  the  use  of  an  electric  motor  system  whose  performance  is  known. 
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The  iterative  procedure  for  defining  14,  is  the  following: 

1.  Provide  an  initial  guess  for  V^. 

2.  Calculate  the  drag  force  of  the  aircraft  and  the  engine  thrust  at  that  value  of  14,.  The  thrust 
of  the  engine  (in  Newtons)  is  given  by 

T  =  -0.0208140  cos  a  +  0.6991,  (5.6) 

where  14  is  in  m/s. 

3.  If  the  component  of  thrust  in  the  drag  direction  is  not  equal  to  the  total  drag  (within  some 
e),  then  define  a  new  14  =  ,  /■■„  Tc°sf^. 

V  aDtotal0.5p5n/ 

4.  Repeat  steps  1  through  3  until  |Tcosa  -  Dtotai|  <  e. 
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5.3  Center  of  Gravity  Analysis 

In  the  next  section,  the  pitching  moment  of  a  given  MAV  design  with  respect  to  the  center 
of  gravity  is  calculated.  Before  that  can  be  done,  however,  the  location  of  the  center  of  gravity 
must  be  defined  first.  The  center  of  gravity  of  the  airplane  is  a  function  of  the  planform,  aspect 
ratio,  area,  and  location  of  the  wing  and  tail  and  the  location  of  the  propeller.  A  schematic  of 
the  procedure  used  to  calculate  the  center  of  gravity  is  shown  in  Figure  5.10. 


Figure  5.10:  Analysis  Procedure:  Center  of  Gravity  Analysis 

The  wing  and  tail  geometry  parameters  are  used  to  calculate  the  weight  of  each  of  these 
surfaces  and  the  location  of  their  centers  of  mass.  Their  position  relative  to  the  origin  is  used  to 
calculate  the  moment  exerted  by  their  weights. 

Also  needed  for  a  mass  balance  analysis  is  the  location  and  weight  of  all  components.  All 
component  locations  are  based  on  the  value  of  the  design  variable  zprop.  Figure  5.11  illustrates 
the  general  layout  of  the  components  in  a  typical  MAV  design. 

The  components  and  placement  constraints  are  dictated  by  the  desire  to  shift  the  center 
of  gravity  as  forward  as  possible,  within  the  constraints  of  the  components  placed  within  the 
fuselage.  The  reason  for  the  desire  of  a  forward  CG  is  that  in  general  the  fuselage  has  to  be 
lengthened  beyond  the  leading  edge  of  the  wing  in  order  to  shift  the  CG  forward  enough  for 


64 
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Figure  5.11:  Dimensions  of  MAV 

a  longitudinally  stable  configuration.  Extending  the  fuselage  increases  the  maximum  dimension 
and  thus  it  should  be  avoided  as  much  as  possible.  Therefore,  by  forcing  the  placement  of  the 
components  as  far  forward  as  possible,  the  required  extension  of  the  fuselage  past  the  leading 
edge  of  the  wing  is  minimized.  The  placement  constraints  of  the  components  are: 

1.  Engine:  The  engine  is  located  just  behind  the  propeller.  A  firewall  (usually  0.3  cm  thick) 
is  placed  directly  behind  the  engine. 

2.  Batteries:  The  two  battery  cells  are  located  just  behind  the  firewall,  leaving  a  small  amount 
of  space  for  foam  cushioning  between  the  firewall  and  the  front  of  the  cells  (about  0.5  cm 
thick). 

3.  Fuel  tank:  The  fuel  tank  is  located  behind  the  firewall,  with  enough  spacing  to  allow  for 
the  fuel  tubing  to  curve  around  through  a  hole  in  the  firewall  towards  the  engine  (about  1 
cm). 

4.  Video  Transmitter:  The  video  transmitter  is  located  at  the  bottom  of  the  fuselage,  be¬ 
hind  either  the  batteries  or  the  fuel  tank  (depending  on  which  of  the  two  extends  further 
downstream).  Foam  cushioning  spacing  is  taken  into  account  as  well. 

5.  Camera:  The  camera  is  located  directly  underneath  the  video  transmitter. 
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6.  Servos:  The  two  servos  are  located  at  the  top  of  the  fuselage,  behind  either  the  batteries 
or  the  fuel  tank  (depending  on  which  of  the  two  extends  further  downstream). 

7.  RC  receiver:  The  radio  control  receiver  is  located  directly  behind  the  servos,  again  allowing 
for  foam  spacing. 

It  should  be  noted  that  the  dimensions,  weight,  center  of  gravity,  and  locational  constraints  of 
each  of  these  components  are  assumed  to  be  known  and  are  furthermore  considered  to  be  constant 
parameters  of  the  design  problem  rather  than  design  variables  (see  Appendix  B  for  details  about 
the  component  specifications).  The  only  true  design  variable  is  rrprop.  All  other  component 
location  variables  are  derived  from  xprop.  Given  this  information  it  is  possible  to  determine  the 
location  of  the  center  of  gravity  of  the  MAY  design. 


5.4  Pitching  Moment  Analysis 

5.4.1  Moment  About  the  CG 

Having  calculated  the  ^-location  of  the  center  of  gravity  of  the  aircraft,  it  is  now  possible  to 
calculate  the  pitching  moment  with  respect  to  the  CG.  The  procedure  is  outlined  in  Figure  5.12. 


Figure  5.12:  Analysis  Procedure:  Pitching  Moment  About  the  CG 
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Referring  to  Figure  5.13,  the  following  expression  for  the  pitching  moment  with  respect  to 
the  center  of  gravity  (positive  for  pitch  up)  can  be  derived: 


-^cg  =  ~(LW  cos  a) (xACw  -  xCg )  -  (Lw  sin  a)(zw) 

-  (LT  cos  a)(xACr  -  xCg )  +  (Ft  sin  d)  (^r) 

-  (Dw  sin  a)(xAcw  ~  ^cg)  +  (Dw  cos a)(zw) 

-  (Dt  sin  a)  (xACt  -xcg)-  ( DT  cos  a)  (zT)  +  (NP) (xprop  +  xCG)  (5.7) 


Figure  5.13:  Pitching  Moment  About  the  Center  of  Gravity 

The  variables  zw  and  zT  represent  the  separation  between  the  centerline  of  the  fuselage  and 
the  quarter-chord  location  of  the  MAC  of  the  wing  and  tail,  respectively.  This  separation  is 
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dictated  by  the  height  of  the  components  that  will  fit  into  the  fuselage  (as  described  in  Appendix 
B).  Furthermore,  as  discussed  in  Section  4.2.3,  the  tail  is  assumed  to  be  always  below  the  wing. 
zw  and  zT  were  both  set  to  2.0  cm  for  use  in  the  optimization  code  of  this  report. 

The  force  NP  is  the  propeller  normal  force  which  is  generated  by  a  rotating  propeller  when 
it  is  inclined  at  an  angle  of  attack  with  respect  to  the  incoming  flow.  Following  the  procedure  of 
Etkin  [23],  the  approximate  value  of  NP  was  calculated  for  the  3-inch  propeller  used  in  the  Cox 
0.010  engine.  At  an  angle  of  incidence  of  30°,  the  normal  propeller  force  was  determined  to  be 
equivalent  to  a  mass  of  1  gram.  In  comparison  with  the  100  gram  takeoff  weight  of  the  MAV, 
this  is  a  small  force  and  it  contributes  only  very  slightly  to  the  pitching  moment.  Therefore  it 
was  decided  to  neglect  the  influence  of  the  propeller  normal  force.  This  simplification  greatly 
shortens  the  computation  time  because  the  pitching  moment  is  no  longer  a  function  of  freestream 
velocity. 


5.4.2  Aerodynamic  Center 

The  location  of  the  aerodynamic  center  of  the  wing  and  tail  also  deserves  some  attention.  As 
discussed  in  Section  3.3,  experiments  of  LAR  wings  at  low  Re  have  revealed  that  the  aerodynamic 
center  of  such  wings  is  not  constant  with  angle  of  attack.  It  generally  shifts  towards  the  trailing 
edge  as  the  angle  of  attack  is  increased.  At  the  time  of  writing  this  report,  the  aerodynamic 
analysis  code  had  not  yet  been  updated  to  include  the  results  of  aerodynamic  center  analysis  of 
Section  3.3.  Future  revisions  of  the  code  will  include  such  analysis  and  will  incorporate  wind 
tunnel  data  into  an  interpolation  scheme  similar  to  the  one  used  in  Section  5.1.1.  As  such,  the 
location  of  the  aerodynamic  center  would  be  a  function  of  angle  of  attack,  aspect  ratio,  and  wing 
planform  and  would  be  determined  empirically  based  on  wind  tunnel  data.  For  the  analysis  code 
used  for  this  work,  a  constant  value  was  assigned  to  the  AC  of  all  wings  and  all  ARs  at  all 
angles  of  attack.  The  value  chosen  was  30%  of  the  MAC. 

For  the  four  wing  planforms  studied  in  this  work,  the  length  of  the  MAC  and  the  ^-location 
of  its  leading  edge  is  calculated  using  the  following  formulas: 

For  rectangular  planforms 


b 


^ root 

MAC 

2/mac 

xmaclb 


Vars 

s 

b 

^root 

0 

0 


(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 
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For  Zimmerman,  inv.  Zimmerman,  and  elliptical  planforms 


(5.13) 

(5.14) 


b  = 

Croot  = 

MAC  = 
Umac  = 

xmacLe  = 


(5.15) 

(5.16) 

(5.17) 


The  geometry  parameters  aa  and  a2  of  the  ellipse-based  wings  are  defined  in  Figure  5.14. 


Figure  5.14:  Geometry  Parameters  for  Ellipse-Based  Wings 
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5.4.3  Longitudinal  Stability  Analysis 

The  longitudinal  stability  requirement  for  the  MAV  design  states  that  the  slope  of  the  pitching 
moment  about  the  center  of  gravity  should  be  negative.  A  test  for  longitudinal  stability  is 
performed  within  the  aerodynamic  analysis  code  by  calculating  the  pitching  moment  at  an  angle 
of  attack  greater  than  the  design  angle  of  attack  of  the  aircraft.  That  is,  the  pitching  moment  of 
the  aircraft  is  taken  at  an  angle  d  +  Act  where  d  is  the  design  angle  of  attack  of  the  aircraft  and 
Ad  is  an  arbitrarily  chosen  increment  in  angle.  For  the  results  presented  in  this  report,  Ad  is  set 
equal  to  2°.  If  the  moment  at  d+ Ad  is  lower  than  the  moment  at  d,  then  the  longitudinal  stability 
requirement  is  satisfied.  The  procedure  used  for  this  stability  analysis  is  shown  schematically  in 
Figure  5.15. 


Interpolate  to  get  x 


Figure 


Analysis  Procedure:  Longitudinal  Stability  Requirement 
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5.5  Maximum  Dimension  and  Penalty  Functions 

The  objective  function  of  a  particular  aircraft  design  can  now  be  calculated.  This  objective 
function  is  a  function  of  the  maximum  dimension  of  the  aircraft  as  well  as  of  the  penalty  functions 
arising  from  the  constraints  of  the  problem.  A  detailed  discussion  of  the  penalty  functions 
associated  with  the  design  problem  is  given  in  Section  4.43.2.  The  steps  followed  for  this  step 
of  the  analysis  code  are  summarized  in  Figure  5.16. 


Maximum  dimension 


Lift  constraint 


Trim  constraint 


Calculate  D„ 


|L-W|<8j? 


Mcg  at  S I  <  e. 


Stall  constraint  (wing)  |  [ct+iw+  ]<  astaI,  w? 


Stall  constraint  (tail) 


Stability  constraint 


Mcg  at  (a+Aa)  <  MCG  at  a 


Objective  Function  f  (x) 


Figure  5.16:  Analysis  Procedure:  Pitching  Moment  About  the  CG 

The  maximum  dimension  of  the  MAV  design  is  calculated  first.  The  distances  in  Figure  5.17 
(dw-y,  dw2,  dri >dr2 >  Hw-t i»  Hw-t2>  Hw,  and  H?)  can  be  determined  based  on  the  geometry  and 
location  of  the  wing  and  tail  and  the  x-location  of  the  propeller.  The  largest  of  these  distances  is 
the  maximum  dimension  of  the  aircraft. 

The  penalty  functions  associated  with  the  constraints  of  lift,  stall,  trim,  and  longitudinal 
stability  are  then  added  to  the  objective  function  after  being  scaled  by  appropriate  penalty  function 
multipliers.  It  is  this  objective  function  value  that  is  to  be  optimized  by  the  genetic  algorithm 
code.  "  ~  ' 
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Figure  5.17:  Dimensions  of  MAY 


5.6  Conclusions  for  Aerodynamic  Analysis  Procedure 

In  this  chapter,  the  procedure  used  to  determine  the  objective  function  of  a  given  set  of  MAV 
design  variables  was  presented.  The  procedure  incorporates  wind  tunnel  data  of  LAR  wings  at 
low  Re  within  an  interpolation  scheme  in  order  to  calculate  lift  and  drag  coefficients  as  a  function 
of  planform,  aspect  ratio,  and  angle  of  incidence.  This  chapter  also  presents  the  procedure  used 
to  determine  the  cruise  speed,  calculate  the  pitching  moment  about  the  center  of  gravity,  and 
perform  a  check  for  longitudinal  stability. 

All  of  these  aerodynamic  performance  characteristics  are  incorporated  into  penalty  functions 
(constraints).  The  maximum  dimension  of  a  given  MAV  design  (a  purely  geometric  performance 
measure)  is  then  combined  with  the  penalty  functions  to  arrive  at  a  final  value  of  the  objective 
function  of  a  particular  design.  This  value  of  the  objective  function  is  then  passed  to  the  genetic 
algorithm  code,  which  attempts  to  optimize  the  design  while  remaining  within  constraints. 
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Chapter  6 

Overall  Conclusions  and  Recommendations 

The  aerodynamic  characteristics  of  wings  of  aspect  ratio  between  0.5  and  2.0  were  experi¬ 
mentally  determined  at  Reynolds  numbers  below  200,000.  Of  primary  interest  to  this  study  were 
the  effects  of  wing  planform  and  aspect  ratio  on  the  forces  and  moments  of  LAR  wings  at  low 
Re.  Results  showed  large  nonlinearities  in  the  lift  curves,  specially  for  aspect  ratios  below  1.25. 
Definite  differences  in  the  performance  of  wings  were  seen  between  the  different  planforms, 
with  the  inverse  Zimmerman  and  rectangular  planforms  being  the  most  efficient  of  the  planforms 
studied.  Analysis  of  the  wind  tunnel  data  included  the  following:  effects  of  Reynolds  number, 
camber,  and  leading  edge  shape,  location  of  the  aerodynamic  center,  and  stall  angle  of  attack  and 
maximum  lift  coefficient. 

The  experimental  data  was  incorporated  within  an  aerodynamic  analysis  procedure  that  used 
both  empirical  and  analytical  techniques  to  determine  the  performance  of  an  aircraft  design.  This 
analysis  procedure  was  made  part  of  a  design-optimization  code  based  on  a  genetic  algorithm  op¬ 
timization  technique.  The  design-optimization  procedure  was  used  to  optimize  an  aircraft  given 
an  example  design  problem.  The  design  problem  required  the  minimization  of  the  maximum 
dimension  of  a  MAV  subject  to  specific  constraints  on  the  lift  required,  trim  conditions,  and 
longitudinal  stability  conditions.  The  genetic  algorithm  code  applied  to  this  problem,  in  conjunc¬ 
tion  with  the  aerodynamic  analysis  procedure  based  on  wind  tunnel  data,  yielded  optimum  MAV 
designs  that  satisfied  all  design  requirements. 

Some  areas  where  further  research  would  be  warranted  have  become  clear  as  a  result  of 
the  present  study.  One  such  area  would  focus  on  interference  effects  of  lifting  surfaces  and 
streamlined  bodies  for  LAR  wings  at  low  Re.  This  particularly  refers  to  a  more  detailed  study 
of  the  effects  of  fuselage  bodies  on  the  drag  and  pitching  moment  of  MAVs.  Also  suggested  is  a 
study  of  the  interaction  of  the  wing  and  horizontal  tail,  effects  of  the  vertical  locations  of  the  two, 
and  effects  of  flow  separation  near  the  trailing  edge  of  wings  and  tails,  specially  in  the  cases  when 
the  control  surfaces  are  located  beyond  the  point  of  separation.  A  more  thorough  understanding 
of  these  aerodynamic  effects  would  greatly  aid  in  the  design  of  more  efficient  MAVs  and  small 
UAVs. 

Another  topic  which  may  be  of  interest  is  a  more  comprehensive  study  of  effects  of  leading 
edge  shape  on  the  lift  and  (perhaps  more  importantly)  pitching  moment  of  LAR  wings.  Prelimi¬ 
nary  results  presented  within  this  report  have  shown  that  the  leading  edge  has  a  dramatic  influence 
on  the  pitching  moment.  Furthermore,  this  influence  is  highly  Reynolds-number-dependent.  A 
study  of  leading  edge  geometries  would  not  only  provide  insight  into  the  most  efficient  leading 
edge  shapes  for  LAR  wings,  but  would  also  be  invaluable  as  a  measure  of  how  much  care  must 
be  taken  in  the  manufacturing  of  MAV  wings. 
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Appendix  A 

Aerodynamic  Results 


A.l  Aerodynamic  Results:  Constant  Aspect  Ratio 


This  section  shows  plots  of  experimental  results  for  constant  aspect  ratio.  That  is,  each  plot 
contains  curves  corresponding  to  each  of  the  four  wing  planforms  at  a  given  aspect  ratio.  These 
graphs  are  useful  to  compare  different  planforms  while  holding  AR  constant.  The  curves  plotted 
are: 


1.  Cl  versus  a  at  Re  =  100,000 

2.  Cd  versus  a  at  Re  =  100, 000 

3.  Cl  versus  Cd  at  Re  =  100, 000 

4.  L/D  versus  a  at  Re  =  100,000 

5.  Cm  versus  a  at  Re  =  140, 000 

6.  hAC  versus  a  at  Re  =  140, 000 
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Figure  Numbers 

Page 

Number 

Description 

AR  =  0.50  | 

A.  1  to  A.4 

78 

Cl  vs.  a,  Cz>  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

AR  =  0.75  | 

A.5  to  A.8 

80 

Cl  vs.  a,  Cd  vs.  a.  Cl  vs.  Cd ,  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

A. 9  and  A.10 

82 

Cm  vs.  a  and  Hac  vs.  ct 

All  four  wing  planforms  at  Re  —  140, 000 

!  AR  =  1.00  | 

A.  11  to  A. 14 

83 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd >  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  —  100, 000 

A.  15  and  A.  16 

85 

Cm  vs.  a  and  Hac  vs.  a 

All  four  wing  planforms  at  Re  =  140, 000 

AR  =  1.25  | 

A.  17  to  A.20 

86 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd >  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

A. 21  and  A.22 

88 

Cm  vs.  a  and  fiAc  vs.  a 

All  four  wing  planforms  at  Re  =  140, 000 

AR  =  1.50  | 

A.23  to  A.26 

89 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

A.27  and  A.28 

91 

Cm  vs.  a  and  Hac  vs.  a 

All  four  wing  planforms  at  Re  ~  140, 000 

AR  =  1.75  | 

A.29  to  A. 32 

92 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd ,  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

A.33  and  A.34 

94 

Cm  vs.  a  and  /iAc  vs.  a 

All  four  wing  planforms  at  Re  =  140, 000 

AR  =  2.00  | 

A. 35  to  A. 38 

95 

Cl  vs.  a,  Cjo  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  four  wing  planforms  at  Re  =  100, 000 

A. 39  and  A.40 

97 

Cm  vs.  a  and  vs.  a 

All  four  wing  planforms  at  Re  =  140, 000 

Table  A.l:  Summary  of  Plots  for  Constant  AR 
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Figure  A.9:  CM  vs.  =  a,  AR  =0.75,  Rec  =140,000 


a  (degrees) 

Figure  A.10:  hA c  vs.  a,  AR  =0.75,  Rec  =140,000 
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a  (degrees) 

Figure  A. 14:  L/D  vs.  a,  AR  =1.00,  Rec  =100,000 
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Figure  A.21:  Cm  vs.  =  a,  AR  =1.25,  Rec  =140,000 


a  (degrees) 

Figure  A.22:  h^c  vs.  a,  AR  =1.25,  Rec  =140,000 
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a  (degrees) 

Figure  A.28:  hxc  vs.  a,  AR  =1.50,  Rec  =140,000 
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a  (degrees) 

Figure  A.29:  Cl  vs.  a,  AR  =1.75,  Re, 


a  (degrees) 

Figure  A.30:  Cp  vs.  a,  AR  =1.75,  Re, 


a  (degrees) 

Figure  A.32:  L/D  vs.  a,  AR  —1.75,  Re, 
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Figure  A.37:  Cl  vs.  Cd,  AR  =2.00,  i?e, 


a  (degrees) 

Figure  A.38:  LID  vs.  a,  AR  =2.00,  Re, 


a  (degrees) 


Figure  A.39:  Cm  vs.  =  a,  AR  =2.00,  Re 


a  (degrees) 

Figure  A.40:  h ac  vs.  a,  AR  =2.00,  Re, 


A.2  Aerodynamic  Results:  Constant  Wing  Planform 

This  section  shows  plots  of  experimental  results  for  constant  planform.  Each  plot  contains 
curves  corresponding  to  each  of  the  seven  aspect  ratios  for  a  given  wing  planform.  These  graphs 
are  useful  to  determine  the  effect  of  varying  aspect  ratio.  The  curves  plotted  are: 

1.  Cl  versus  a  at  Re  =  100,000 

2.  Cd  versus  a  at  Re  —  100, 000 

3.  Cl  versus  Cd  at  Re  =  100, 000 

4.  L/D  versus  a  at  Re  =  100,000 

5.  Cm  versus  a  at  Re  =  140, 000 

6.  h\c  versus  a  at  Re  =  140, 000 


Figure  Numbers 

Page 

Number 

Description 

Rectangular  Planform 

A.41  to  A .44 

99 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  seven  aspect  ratios  (0.50  <  AR  <  2.00)  at  Re  =  100, 000 

A.45  and  A.46 

101 

Cm  vs.  a  and  Hac  vs.  a 

Six  aspect  ratios  (0.75  <  AR  <  2,00)  at  Re  =  140, 000 

Zimmerman  Planform 

A. 47  to  A.50 

102 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd ,  and  L/D  vs.  a 

All  seven  aspect  ratios  (0.50  <  AR  <  2.00)  at  Re  =  100, 000 

A.51  and  A.52 

104 

Cm  vs.  a  and  /iac  vs,  a 

Six  aspect  ratios  (0.75  <  AR  <  2.00)  at  Re  =  140, 000 

Inverse  Zimmerman  Planform 

A.53  to  A.56 

105 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  seven  aspect  ratios  (0.50  <  AR  <  2.00)  at  Re  =  100, 000 

A.57  and  A.58 

107 

Cm  vs.  a  and  /iac  vs.  a 

Six  aspect  ratios  (0.75  <  AR  <  2.00)  at  Re  =  140, 000 

Elliptical  Planform 

A.59  to  A.62 

108 

Cl  vs.  a,  Cd  vs.  a,  Cl  vs.  Cd,  and  L/D  vs.  a 

All  seven  aspect  ratios  (0.50  <  AR  <  2.00)  at  Re  =  100, 000 

A.63  and  A.64 

110 

Cm  vs.  a  and  h ac  vs.  a 

Six  aspect  ratios  (0.75  <  AR  <  2.00)  at  Re  =  140, 000 

Table  A.2:  Summary  of  Plots  for  Constant  Wing  Planform 
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a  (degrees) 

Figure  A.45:  Cm  vs.  =  a,  rectangular  planform,  Rec  =140,000 


a  (degrees) 

Figure  A.46:  /iac  vs.  a,  rectangular  planform,  Rec  =1 


EH 


a  (degrees) 

Figure  A.51:  Cm  vs.  =  a,  Zimmerman  planform,  Rec  =140,000 


a  (degrees) 

Figure  A.52:  h ac  vs.  a,  Zimmerman  planform,  Rec  =140,000 
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8 


a  (degrees) 

Figure  A.63:  Cm  vs.  =  a,  elliptical  planform,  Rec  =140,000 


a  (degrees) 

Figure  A. 64:  h^c  vs.  a,  elliptical  planform,  Rec  =140,000 
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Appendix  B 
Components 


This  section  discusses  the  components  chosen  for  the  MAV  design  problem  considered  in 
this  report.  The  components  listed  here  are  the  smallest  and  lightest  off-the-shelf  components 
available  at  the  time  of  this  writing.  In  the  dimension  specifications  that  follow,  the  “length”  is 
taken  to  be  in  the  x  direction,  the  “width”  is  taken  to  be  in  the  y  direction,  and  the  “height”  is 
taken  to  be  in  the  z  direction. 

1.  Engine  and  Propeller:  The  engine  chosen  for  the  MAV  is  a  Cox  glow  engine  with  a 
displacement  of  0.164  cubic  centimeters  (0.010  cubic  inches).  It  uses  a  7.62  cm  (3  in) 
propeller.  It  runs  at  full  throttle  for  the  entire  flight  as  it  does  not  have  a  carburetor. 


Dimensions  3.0  cm  long  x  1.5  cm  wide  x  3.5  cm  high 

Mass  13.6  grams 

Center  of  Gravity  (0.50  x  length)  cm  from  front  edge 


2.  Fuel  Tank:  The  fuel  tank  can  hold  a  maximum  of  6  cubic  centimeters  of  fuel  (9  grams  of 
fuel).  This  allows  the  engine  to  run  for  approximately  10  minutes. 


Dimensions 

Mass 

Center  of  Gravity 


4.3  cm  long  x  1.7  cm  wide  x  1.7  cm  high 
11.0  grams 

(0.50  x  length)  cm  from  front  edge 


3.  Video  Transmitter:  The  video  transmitter  chosen  is  a  2.4  GHz  FM  transmitter  with  an 
output  power  of  80  mW.  It  is  capable  of  transmitting  a  color  video  image  up  to  two  miles 
away. 


Dimensions  5.5  cm  long  x  2.0  cm  wide  x  0.75  cm  high 

Mass  7.4  grams 

Center  of  Gravity  (0.50  x  length)  cm  from  front  edge 
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4.  Video  Camera:  The  video  camera  is  a  color  CMOS  camera  with  a  pinhole  lens. 


Dimensions  1.5  cm  long  x  2.3  cm  wide  x  1.5  cm  high 

Mass  7.0  grams 

Center  of  Gravity  (0.50 x  length)  cm  from  front  edge 


5.  Radio  Control  Receiver:  The  RC  receiver  is  a  micro  6-channel  receiver  (FMA  Extreme). 


Dimensions  4.7  cm  long  x  2.3  cm  wide  x  1.5  cm  high 

Mass  10.6  grams 

Center  of  Gravity  (0.50  x  length)  cm  from  front  edge 


6.  Radio  Control  Servos:  The  servos  are  sub-micro  servos  (Hitec  HS-50). 


Dimensions  2.0  cm  long  x  1.2  cm  wide  x  2.5  cm  high 

Mass  5.6  grams 

Center  of  Gravity  (0. 50 x  length)  cm  from  front  edge 


7.  Batteries:  The  batteries  power  both  the  video  equipment  and  the  radio  control  equipment. 
Two  lithium  cells  are  used  wired  in  series.  They  have  a  capacity  of  300  mAh  and  could 
run  the  electronics  in  the  MAV  for  a  duration  of  approximately  30  minutes. 

Dimensions  3.0  cm  long  x  0.5  cm  wide  x  2.5  cm  high 

Mass  7.5  grams  each 

Center  of  Gravity  (0.50x length)  cm  from  front  edge 


112 


